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Abstract: 

We study the scale dependence of the twist-3 quark-gluon parton distributions using 
the observation that in the multi-color limit the corresponding QCD evolution equations 
possess an additional integral of motion and turn out to be effectively equivalent to the 
Schrodinger equation for integrable open Heisenberg spin chain model. We identify the 
integral of motion of the spin chain as a new quantum number that separates different 
components of the twist-3 parton distributions. Each component evolves independently 
and its scale dependence is governed by anomalous dimension given by the energy of the 
spin magnet. To find the spectrum of the QCD induced open Heisenberg spin magnet 
we develop the Bethe Ansatz technique based on the Baxter equation. The solutions 
to the Baxter equation are constructed using different asymptotic methods and their 
properties are studied in detail. We demonstrate that the obtained solutions provide a 
good qualitative description of the spectrum of the anomalous dimensions and reveal a 
number of interesting properties. We show that the few lowest anomalous dimensions 
are separated from the rest of the spectrum by a finite mass gap and estimate its value. 
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1. Introduction 



The evolution equations play an important role in the QCD studies of hard processes as 
they allow to find the dependence of hadronic observables on the underlying high-energy 
scales It has been recognized recently that in different kinematical limits the QCD 
evolution equations possess an additional hidden symmetry. Namely, the Regge asymp- 
totics of hadronic scattering amplitudes H and the scale dependence of the leading 
twist light-cone baryonic distribution amplitudes 0] reveal remarkable properties of inte- 
grability. Both problems turn out to be intrinsically equivalent to the Heisenberg closed 
spin magnet which is known to be integrable (1 -|- 1)— dimensional quantum mechanical 
system and their solutions can be found by applying powerful methods of Integrable 
Models [§. In the present paper we continue the study of integrable QCD evolution 
equations initiated in |^ by considering the evolution equations for twist-three quark- 
gluon distribution functions. These functions describe the correlations between quarks 
and gluons in hadrons and provide an important information about the structure of 
hadronic states in QCD. They determine the twist-3 nucleon parton distributions [0 
and the high Fock components of the meson wave functions 0] which in turn can be 
accessed experimentally through the measurement of different asymmetries. 

The twist-3 quark-gluon distribution functions D{xi, X2, X3] /i) can be defined in QCD 
in terms of hadronic matrix elements of nonlocal gauge invariant light-cone operators 
J^{zi,Z2,Z3) as HQ 

{p,s\J^{z,,Z2,z,)\p,s) = ^'vxe'^^^'^+^'''+^^'^^^P^^D{x^,X2,x,■fi), (1.1) 



where Vx = dxidx2dx2,5{xi + X2 + x^) and |p, s) is the nucleon state with momentum 
p and spin s. The distribution function D{xi^X2,x^; ^i) vanishes outside the region 
— 1 < a;i^2,3 < 1 and the scaling variables Xi, X2 and — X3 have a meaning of hadron 
momentum fractions carried by quark, gluon and antiquark, respectively, in the infinite 
momentum frame 0. For < a;j < 1 (or — 1 < a;j < 0) the corresponding parton belongs 
to the initial (or final) state nucleon. The matrix element entering the l.h.s. of (|1.1| ) 
contains ultraviolet divergences that are subtracted at the scale /x. The normalization 
scale /i defines the transverse separation of the partons and is determined by the hard 
scale of the underlying process [|l|. 

The operator J^{zi, Z2, Z'i) describes the interacting system of quark, antiquark and 
gluon "living" on the light-front = zn^ with = and separated along the 
z— direction. The twist-3 operator J^{zi, Z21 z^) is given by one of the following expres- 



sions ITDI, ITT 



S ^{zi, Z2, z^i) = q{zin) G ±y{z2n) ± iG ^^{z2n)'^^ n^flqizsn) (1.2) 
Tr{zi,Z2,Z3) = q{zin)n^a''fn''G^p{z2n)T qiz^n) (1.3) 



^The twist-3 contribution to the wave function is given by a similar expression with xi + X2 + = 1 
and the matrix element calculated between the vacuum and the meson state. 
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with r = {11,275}, G^u = being dual gluon field strength and "±" denoting 

the "transverse" Lorentz components orthogonal to the plane defined by the vectors 
Pfj^ and nfj_. The definition of the chiral-even and chiral-odd quark-gluons distributions 
corresponding to the operators ( [1.2|) and ( |1.3|) , respectively, and their relation to the 
twist-3 nucleon parton distributions are given in the Appendix A. In order to avoid 
additional complication due to mixing of the chiral-even operators with pure gluonic 
operators we will assume quarks to be of different (massless) flavor. It is implied that the 
gauge invariance of the operators ( |1.2|) and (|1.3| ) is restored by including non-Abelian 
phase factors Pq"^! '^^''^t^i^) connecting gluon strength fleld with quark and antiquark 
flelds. 

In this paper we shall study the /x— dependence of the distribution functions ( |1.1D in 
the multi-color QCD. The standard way of flnding this dependence O IT3 consists 



in expanding nonlocal operator J^{zi, Z2, Z3) entering ( |1.1| ) over the set of local gauge 

invariant composite operators ON,n = Qin D)^G{n D)^~'^q with n = 0, + 1 and 

D^j, = ~ ig^ti being the covariant derivative. In this way, one flnds the evolution of 
the moments J^^dxidx2dxs5{xi + X2 + x-^) Xi~"'x'^ D(xi) by diagonalizing a nontrivial 
(A^ -|- 1) X (A^ + 1) mixing matrix of the operators O^^n- It becomes straightforward to 



calculate numerically the spectrum of the anomalous dimensions for any given A^ ||T3 
but the general analytical structure of the spectrum remains unknown. To overcome 
this problem we apply the approach developed in ^j. It allows to construct the basis 
of local twist-3 quark-gluon operators which evolve independently and whose anomalous 
dimensions can be calculated analytically. 

According to the deflnition, Eq. (p..l|), the /i— dependence of the twist-3 quark-gluon 
distributions follows from the renormalization group evolution of the corresponding non- 
local operators (|1.2|) and ([L.3|). To leading logarithmic accuracy one flnds illl 



^ -J^{Zi, Z2, Z3) = -^^^Hjr ■ J^{zi, Z2, Z3) (1.4) 



din fx 27r 

with TijF being some integral operator describing a pair- wise interaction between quark, 
gluon and antiquark with the light-cone coordinates Zi, Z2 and z^, respectively. The 
approach of is based on the identiflcation of the evolution equation ( |1.4| ) as one- 
dimensional Schrodinger equation for three particles with a pair-wise interaction. The 
light-cone coordinates of the quarks and gluon = ZiU^ deflne the spacial coordinates, 
Zi, whereas the evolution parameter In/x^ plays the role of the evolution time. The inte- 
gration of three-body Schrodinger equation (|1.4|) becomes problematic unless there exists 
an additional symmetry. It turns out that this is the case for the twist-3 quark-gluon 
correlation functions D{xi;fi). The corresponding evolution equation ( |1.4| ) possesses an 
additional "hidden" integral of motion in the multi-color limit and, as a consequence, 
it is completely integrable [Q. As we will show, the underlying integrable structure is 
identical to that for the integrable open noncompact SL{2] M) spin chain models. Sim- 
ilar conclusions have been reached in the recent publications [|14|. Our results partially 



overlap with the results of and we shall comment on their relation below. 

"Noncompactness" of the spin group is a novel feature that QCD brings into the 
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theory of integrable open spin magnets and that makes the analysis of the evolution 
equations interesting on its own. We would like to notice that an additional motivation 
for studying noncompact integrable open spin chains comes from the recent analysis of 
the Regge asymptotics of quark-gluon scattering amplitudes We also expect that 



these models (in their generalized version) will inevitably appear in the analysis of the 
evolution equations for higher twist quark-gluon correlation functions 0. 

In the present paper we develop an approach for solving noncompact open spin chains 
models based on the Baxter equation and the fusion relations and apply it to find the 
spectrum of the anomalous dimensions of the twist-3 quark-gluon correlation functions in 
the leading logarithmic approximation in the multi-color QCD. As we will demonstrate 
this spectrum exhibits a number of interesting properties. The few lowest anomalous 
dimensions turn out to be separated from the rest of the spectrum by a finite "mass gap" 
and can be calculated exactly. The remaining part of the spectrum can be calculated 
using different asymptotic methods and the obtained approximate expression agree well 
with the exact results for large spin A^. 

The presentation is organized as follows. In Section 2 we explore the conformal 
symmetry of the evolution equations to write the QCD evolution kernels in the "normal" 
form in which their conformal invariance becomes manifest. In Section 3 we construct 
a general integrable open spin magnet and identify the values of parameters for which 
the Hamiltonian of the model coincides with the QCD evolution kernels for nonlocal 
operators and T. Section 4 is devoted to diagonalization of the conserved charge by 
means of the Algebraic Bethe Ansatz. Section 5 contains a derivation of the expression 
for the energy of the open spin magnet, or equivalently the spectrum of the anomalous 
dimensions of the quark-gluon correlation functions, in terms of the solutions to the 
Baxter equation. In Section 6 we study the properties of the Baxter equation and 
find its exact solutions. In Section 7 the asymptotic methods are applied to describe 
a fine structure of the spectrum. Summary of the main results is given in Section 8. 
In Appendix A the relation between the quark-gluon distributions and twist-3 nucleon 
structure functions is discussed. Appendix B contains definition and some properties of 
the Wilson polynomials. 



2. QCD evolution kernels 

In the multi-color limit the QCD evolution kernel Tijr entering ( |1.4| ) gets contribution 



only from planar diagrams and to one-loop level it is equal to the sum of the quark- 
gluon and antiquark-gluon interaction kernels while the interaction between quark and 
antiquark is suppressed by a color factor l/N^. The kernel TCj^ depends on the choice of 
nonlocal operator and for twist-3 nonlocal operators defined in ( |1.2|) and ( |1.3|) it is given 
in the multi-color limit by the following expressions ^ 

^T = ^(j,, + +i'{j,9-l) +^{j9^+f) +^ (2-1) 
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with = d\nT{x)/dx, 6 = — 3/2 + and Jik being the operators acting on the 
hght-cone coordinates of the partons {i,k = q,g,q) and defined below in ( |2.5|) . Here, 
the kernel TYt governs the RG evolution of the operators Tr(-2j). To one-loop level it 
does not depend on the choice of T and we will not indicate the subscript T explicitly. 
The kernel Ti-s- is obtained from T-Cs+ by interchanging the coordinates of quark and 
antiquark and therefore it will not be considered separately. 

It is well known that to one-loop level the QCD evolution kernels ( |2.1| ) and 



( p.2|) inherit the conformal symmetry of the bare QCD Lagrangian which is reduced 
on the light-cone to its SL{2,M.) subgroup generating a projective transformation z — > 
{az + b)/{cz + d) on the line = zn^ with ad — be = 1. Denoting the SL{2,M.) 
generators as L+, L_ and Lq we define their action on the quark and gluon fields ^{z) = 
{q{nz), q{nz), n^G^y{nz)} on the light-cone = as 

L_,fc$(2;fc) = -d,^^{zk) 

L+^k^Zk) = {zld,, + 2jkZk) Hzk) (2.3) 

Lo,fc$(zfc) = {zkdz^^ +jk) ^(zk) , 

where jk = {h + Sk)/2 is the conformal spin of the field $fe {k = q,q,g) defined as 
the sum of its canonical dimension, 4, and the projection of its spin on the line zn^, 
= isk^k- In particular, using lq = Iq = 3/2, Ig = 2 together with Sg = Sg = 1/2, 
= 1 we get the conformal spins as 

jg = J, = 1 , jg = 3/2 (2.4) 

for (anti-) quark and gluon fields, respectively. 

The operators Jqg and Jqg entering ( |2.2D and ( |2.1| ) are defined in terms of two-particle 
Casimir operators as follows 

^ik ~ Jik{Jik ~ 1) ; = + -^fc.a (2-5) 

with i^k = q ,g ,q and a = 0, ±. The eigenvalues of the operator Jqg give the possible 
values of the spin in the quark-gluon channel and have the form jqg = jq+jg+n = 5/2+n 
with n being a nonnegative integer. The operator Jgg has a similar interpretation in the 
gluon- antiquark channel. The SL{2) symmetry of the evolution kernels (|2.1|) and (|2.2| ) 
becomes manifest 

[Ut, L„] = [Ut, L'] = [L^, L']=0, (2.6) 
where La (a = ±, 0) are the total three-particle SL{2) generators 

La = La^q + + ; ~ 2 ~^ ^-^+) "^0 • (2-7) 

Note that the conformal invariance holds separately for all terms entering the 1/iVc— expan- 
sion of the evolution kernels (12. 21) and (E.ll). 
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Going from the matrix element of nonlocal operator to the quark-gluon distribution 
function, Eq. ( |1.1| ), we find that D{xi,X2,Xs) obeys the RG equation similar to ( |1.4| ) 
with the Hamiltonian Tijr transformed from the coordinate z— representation into the 
momentum x— representation. The analysis of the evolution equation for D{xi,X2, x^) is 
based on the solutions to the Schrodinger equation 0] 

^-J^ ■ '^N,q{Xi,X2, X3) = SN,q'^N,q{^l, X2, X3) , (2.8) 

where the eigenstates \E'jv ,j(a;j) are homogeneous polynomials in Xi of degree N and index 
q enumerates different energy levels SN,q- Here, the Hamiltonian is given by the evolution 
kernels and (|2.2| ) with the S'L(2;M) generators defined in the x— representation as 

Lkfi'^ixi) = {xkdk+ jk)^ixk) 

Lk,+^{xi) = (2.9) 
Lfc,„^(xi) = {xkdl^ + 2jkd^^)^{xk) . 



Finally, solving the Schrodinger equation ( p.8| ) one can construct the basis of local con- 
formal operators 

Civ,g(0) = ^A^,g(^9^l,^9^2,^9^3)J^(-^l,-^2,-^3) , (2.10) 

Zi=0 

which have a fixed operator dimension, or equivalently a fixed total number of covari- 
ant derivatives, and which do not mix under renormalization. It follows from (|2.8|) that 
their anomalous dimensions are determined by the energy of the state '^N,q 

7.(iV) = ^^N,q (2.11) 

while their matrix elements over nucleon states evaluated at a low normalization scale 
fiQ define the set of nonperturbative dimensionless parameters 

07V,g(/^o) = {P, s\ON,g{0)\p, s) ^=^^ . (2.12) 

Substituting (|1.1|) into ( |2.1CI| ) one finds the solution to the evolution equation for (gen- 
eralized) moments of the distribution functions as 



Vx-^N,q{Xl,X2,X3)D{xi,X2,X3) = (f)N,q{lJ'o) [ 7 T , (2.13) 

where 60 = ll^c/3 — 2nj/3 and Vx = dxidx2dx3S{xi + X2 + X3). Inverting this relation 
one can obtain the expression for the distribution function as an integral over complex 
moments A^. 

■^In the case of the twist-3 quark-gluon wave function, xi + X2 + = 1 and Xi > 0, the relation 
(2.13) can be inverted using the completeness condition for the set the states 'i'N,q{xi) with N — 
nonnegative integer as 

D{xi,X2,X3; fi) = Xixlxs 'l)N,q'^ N.q{xi) I —. r I 



N>0,q 
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Thus, the scale dependence of the twist-3 quark gluon distributions becomes effec- 
tively equivalent to the eigenproblem ( |2.8| ). In what follows we shall solve exactly the 
Schrodinger equation (|2.8|) and study in detail the properties of its solutions. Analyzing 
( p.8|) one has to identify the total set of the conserved charges. The conformal symme- 
try allows to identify the degree N of the polynomials '^N,q{xi) as a trivial integral of 
motion. For fixed the eigenstates \E'Ar,g(xj) belong to the irreducible representation of 
the SL{2, M) group parameterized by a total conformal spin h = N + jg + jg + jg and 
satisfy the relations 

7 

Lo'^hixi) = h^h{xi) , h = N + - 

L_^hixi) = (2.14) 
L^^hix^) = h{h-l)^h{xi). 

Thus defined "ifhixi) can be identified as the highest weights of the SL{2,'R) represen- 
tations labeled by a nonnegative integer N. Note that for given "^hixi) there exists an 
infinite series of the "excited" eigenstates, (L+)™\E'/i = (— xi —X2 —X3)"^'ifh, which vanish 
due to Xi + X2 + Xs = 0. 

However the conformal symmetry does not fix the eigenfunctions uniquely. 
This can be seen in a number of ways. According to ( p. 71) the total conformal spin h 
is equal to the sum of the conformal spin of the quark-gluon system and the conformal 
spin of antiquark. For given h the value of the spin in the g^f— channel can be arbitrary, 
jqg = jq + jg + = 5/2 + 171 with < m < A^, and this leads to an additional degeneracy 
in solutions to ( 2.14| ). The degeneracy can be removed if there exists one more SL{2) 



invariant conserved charge Q. In this case, the Schrodinger equation (|2.8|) becomes 
completely integrable - the number of degrees of freedom matches the number of the 
conserved charges, Lq, and Q, and, as a consequence, the wave functions, '^N,q, and 
the energy, SN,q, are uniquely fixed by the conformal spin A^ and the eigenvalues, q, of 
the charge Q. 

It turns out that such operator exists for the evolution kernels corresponding to the 
twist-3 quark-gluon operators in the multi-color limit, Eqs. ( p.l|) and (|2.2| ) 



[ns±,Qs^] = o, [nT,QT] = o. (2.15) 

The SL{2) invariant charges Qs+ and Qt are given by 

Qs^ = {Ll,,Ll,}-hl,-hl, (2.16) 

2 2 9 2 9 2 

Qt = {Li2, L23} - ~ 2^23 y (2-17) 

where { , } stands for an anticommutator, = L^g and L23 = L'^g are two-particle 
Casimir operators in the quark-gluon and gluon-antiquark channels. The charge Qs- is 
obtained from Qs+ by permutation of the particles (123) (321). The kernel TYt is 
invariant under this transformation and, as a consequence, its eigenstates have a definite 
parity under zi ^ ,23. 
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3. Integrability of the effective QCD Hamiltonian 



Let us demonstrate that the Schrodinger equation (|2.8| ) defined by the evolution kernels 
Eqs. (|2.1|) and (|2.2| ) is equivalent in the multi-color limit to the one-dimensional Heisen- 
berg open spin chain model. To this end, we shall apply the Quantum Inverse Scattering 
Method m, to construct an integrable (inhomogeneous) open spin chain model and 



show that for certain values of parameters the Hamiltonian of this model coincides with 
the QCD evolution kernels (O) and (O). 



3.1. Open spin chains 

To start with we consider M particles with the coordinates x„ (n = 1, ...,M) on a line 
and assign to each particle three SL{2) generators (i^+,n5 -^-,n, -^o,n)- In what follows we 
shall refer to these generators as to spin— j„ operators. They are realized as differential 
operators (|2.3|) acting on the quantum space of the n—th. particle that we denote as V^„. 
The interaction between M particles occurs through a pair-wise interaction between 
the nearest-neighbor spins. The corresponding pair-wise Hamiltonians depend on the 
two-particle Casimir operators 

5^ (^aj + , J,A; = 1,...,M (3.1) 

Q=0,l,2 

with L-tj = Li j ± zL2j. 

The definition of the model is based on the existence of the solution to the Yang- 
Baxter equation 

Rhh iu)Rj,j, {v)Rj^j, {u~v) = Rj^j^ {u - v)Rj,j, {v)Rj,j., (u) . (3.2) 



Here, the operator Rj^j2 (u) acts on the tensor product of the quantum spaces Vj 



32 



and depends on an arbitrary complex spectral parameter u. The solution to ( |3.2| ) is 
given by 0, 

Rnn(u) = (-)™ ''t':'''!''^'t^'^'"' . (3.3) 

r(Ji2 + lu) r(ji + j2 - 

where the operator J12 is given by the sum of two SL{2) spins and it is defined as a 
formal solution to 

Jl2(Jl2-l) =i^?2- (3.4) 

For arbitrary SL{2) spins ji and j2 the eigenvalues of the operator J12 have the form 
J12 = ji + j2 + n with n = 0,1,..., which correspond to decomposition of the tensor 
product Vj^ over the irreducible SL{2) components of spin ji2- 

Having the explicit expression for the i?— operator, Eq.( p.3|) , we construct an in- 
tegrable open spin model following Sklyanin Namely, we define the monodromy 
operator acting on the space Vj ®Vj^ ® ... ® Vj^j 

Tj{u) = Rj^j{u — iuJi)Rj2j{u — iuj2)...Rjjyjj{u — iujm) (3.5) 
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with j being the spin of the auxihary space Vj and Uk being the shifts (spin chain 
impurities) associated with the A;— th particle. Let us impose the additional condition 
that the spins and the impurities of all "intermediate" sites are the same 

J2 = J3 = ••• = JM-I , Ll!2 = UJ3 = ... = UJm-I , (3.6) 

while j2 7^ JiJm and uj2 7^ uji,ujm- 

As we will show in Sect. 3.2 it is this condition (together with 002 = and M = 3) 
that one finds matching the QCD evolution kernels, ( p.l|) and ( p.2|) , into integrable spin 
chain Hamiltonian. Then, the transfer matrix of the open spin chain with M sites is 
defined as 

t,{u)=ti,[T,{u)Tr\-u)] . (3.7) 

Here, the trace is taken over the auxiliary SL{2) representation space of the spin j. 
Choosing different values of j one obtains (an infinite) set of the transfer matrices tj{u) 
which by virtue of the Yang-Baxter equation (|3.2|) commute with each other for different 
values of the spectral parameters 

[t,,{u),t,,{v)] = (3.8) 

for any ji and j2- As a consequence, expanding tj{u) in powers of the spectral param- 
eter u and choosing different values of the spin j one obtains the family of mutually 
commuting operators. Their explicit form depends on the set of parameters, jk and Uk, 
which define the spectrum of spins and impurities of the model, respectively. This family 
contains the Hamiltonian 7i^/ of the model as well as a complete set of the conserved 
charges. In addition, the transfer matrix possesses an additional SL{2) symmetry]^ 

[tj{u),La] = , La = La,i + ... + La,M , a = 0,±, (3.9) 
which ensures the SL{2) invariance of the Hamiltonian, (|2.6| ). 



3.2. Integrals of motion 



Let us show that the QCD integrals of motion ( p.l7|) and ( |2.16| ) appear in the expansion 
of the spin— 1/2 transfer matrix, t_i/2{u), given by ( |3.7| ) with the auxiliary SL{2) spinQ 
j = —1/2 and the number of sites M = 3. 

%, that the spin— 1/2 i?— operator (p. 3]) coincides with the Lax 



It is well known 



operator of the Heisenberg spin chain, R_i 



Lj (u) 



u -\ i 

2 



=o,± 



u + {-iL^ 



-iL^ 



u 



-iL- 



(3.10) 



"'General definition of integrable open spin chain models involves the boundary matrices K±{u) 



We have chosen the simplest solution K± 



satisfying the reflection Yang-Baxter equation |2^, ^ | 
to ensure the SL{2) invariance of the Hamiltonian, Eq.(2.6|). 

^Indeed, it follows from the definition (3.5) that the monodromy operator commutes with the sum 
of the total spin L and the auxiliary spin Lj, \Tj{u), Lj + L] ~ 0. To get (^) one has to take the trace 
over Vj and use its cyclic symmetry property. 

^Throughout this paper we shall use the convention jsu(2) ~ ~is 



J SL{2) ■ 
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with La being spin— j generators and L± = Li ±iL2. Substituting ( |3.10| ) into ( pT5| ) one 
evaluates the spin— 1/2 monodromy operator as a 2 x 2 matrix 

T_i/2(m) = - zcJi)...Lj,,(u - iuM) = (^^1^1 ^1^1 ^ . (3.11) 

Here, a{u), b{u), c{u) and d{u) are expressed in terms of the spin operators and satisfy 
the Yang-Baxter relations of the form 

[6(«),c(^;)] = ^^M^M^^(!M!^) (3.12) 

together with similar relations for other components of T_i/2(m)- The spin— 1/2 transfer 
matrix t^i/2{u) is given by (p.7|) with T_i/2(m) replaced by its expression ( |3.11| ). It 
proves convenient to change the normalization of the transfer matrix by introducing the 
following operator [|17 



t-i/2iu) = t_i/2y-u--jp{u), (3.13) 
where the c-valued factor p{u) compensates the poles of t_i/2(— u — and is given by 

M 

p{u) = J](n + ij^){u - ZJ+ + i) (3.14) 

k=l 

with the parameters defined as 

jk=jk±i^k- (3.15) 



Thus defined normalized transfer matrix can be expressed as [|17 



f_l/2(M) = trfLi/2(M) , f_i/2(M) = r_i/2|^M - (T2Tli/2(^-^i - Cr2 (3.16) 

and it has the following properties. Firstly, in contrast with t_i/2{u), the normalized 
transfer matrix does not have poles in u and, secondly, t-i/2{u) has the form of a poly- 
nomial of degree M in u"^ with the corresponding operator coefficients providing the set 
of the conserved charges of the model. 

It becomes straightforward to calculate the normalized transfer matrix t_i/2(^i) for 
the chain with M = 3 sites. Substituting (|3.1CI|) into (|3.11|) and (|3.16|) one finds after 
some algebraQ 

tli/2(n) = (W + l)r(n2)-2 J] {u'~LI + luI) (3.17) 



fc=l,2,3 



^Similar calculation of the transfer matrix has been performed in the second publication in p4[ |. We 
disagree with Eq. (47) there in the coefficient in front of r(u^). 
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with one-particle Casimir operators Lj, = jk{jk — 1) being c- numbers and the operator 
r(M^) defined as 

riu')=u'[L'- LI]+L^[lI + uI]-\q 

fc=l,2,3 

+LlLl - Ljooj - Lliujf - ooi + ool) - Ljujf . (3.18) 

The charge Q is given by 

Q = -2^2 ^23]- + {Ll2, ^23}+ + 2(^2' - ^3) + 2iu;l - cuf) L^, . (3.19) 

with L"^!^ and being two- and three-particle Casimir operators defined in ( p.l|) and 
( p.7|) , respectively. 

The normalized transfer matrix (|3.17|) involves two mutually commuting operators, 
Q and L^. Together with Lq they form the complete set of the conserved charges for 
integrable open chain of M = 3 particles with arbitrary values of the spins jk and the 
impurity parameters Uk- To fix their values we compare ( |3.19| ) with the expressions for 
the integrals of motion of the QCD evolution equations, Eqs. ( p.l7|) and ( p.l6|) . We find 
that the spins of particles are equal to the conformal spins of the corresponding fields 

ji = is = 1 , J2 = ^ , (3.20) 
and the impurity parameters are given by 

Qs+: ^? = |> ^2 = 0, ul = ^ (3.21) 

Qt : ^1 = ^, ^2 = 0, ujI = ^. (3.22) 

These matching conditions define the parameters Ui and up to a sign. As we will 
show in Sect. 3.3, the Hamiltonian of the model is invariant under cji 3 — c^i,3. 

Eqs. ( |3.20|) , ( p.21| ) and (|3.22| ) establish the equivalence relations between the QCD 



and 



evolution kernels and the open spin chain models. Instead of considering two cases 
( p.21| ) and ( p.22| ) separately, we shall treat ui and as free parameters and introduce 
the charge 

Q(^i, ^3) = {Ll„ Ll,}+ - 2ujI Ll, - 2ul Ll, . (3.23) 

In what follows we shall find the fiow of the eigenvalues of this charge, q{uJi,uj3), in the 
parameters ui^s and identify the integrals of motion, qt and qs+, as corresponding to the 
special values of the fiow parameters 

,,_,(|i), ,.^,(y). (3.24) 
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3.3. Integrable Hamiltonian 

Let us complete an identification of tlie QCD evolution equations as integrable systems 



by showing that the evolution kernels ( |2.1|) and ( |2.2| ) coincide with the Hamiltonian of 



open spin chain model. The latter enters into the expansion of the transfer matrix tj 



u 



given by ( p.7|) with the auxiliary spin j equal to the spin of intermediate "gluonic" site 
j = j2 = 3/2. The corresponding monodromy operator takes the form 



Tj{u) = Rj,j{u - iui)Rj^j{u)Rj,j{u - iu^] 



(3.25) 



3=32 



where uj2 = 0, ji = js = 1 and j2 = 3/2. 

We would like to stress that the spin chain described by ( ^.25| ) is inhomogeneous - the 
spins and impurities of the end-points are different from the ones of the internal sites of 
the chain. To our best knowledge such QCD induced spin magnets have not been studied 
before and their analysis represents a certain interest from point of view of integrable 
models. It is for this reason that we generalize the definition ( p. 25 ) by enlarging the size 
of the open spin chain to an arbitrary number of sites, M > 3, and fixing the parameters 
(spins and impurities) according to 



Tj{u) = Rj,j{u 



where we put U2 



iuJi)Rj,j (u) . ..Rj^j_^j {u) Rj,,,j {u 



(3.26) 



3=32= 



=3M-l 



ujm-1 = in order to match (|3.25| ) for M = 3 



Finally, one finds the transfer matrix tj{u) by substituting the monodromy operator 
( p.26| ) into ( |3.7| ) and defines the Hamiltonian of inhomogeneous spin chain with M sites 



as 



n 



M 



-dy,lntj{u 



0) = -9Jntr [T^{u)Tr\-u)] 



u=0 



(3.27) 



We would like to notice that this relation differs from the standard definition of Hamilto- 
nian based on the so-called fundamental monodromy operator [0, |22|. In contrast with 
the latter, the spins of quantum spaces Vj-^ and Vjj^^ corresponding to the end-points of 
the chain ( p.26| ) are different from that of the auxiliary space Vj. f\ 

One finds the explicit form of the Hamiltonian ( |3.27| ) in the standard way [|T^] by 
replacing the monodromy operator Tj (u) by its expression ( |3.26| ) and taking into account 
the following property of the i?— operator ( pl3| ) 



%.(0)=P, 



3k3 ' 



k 



M-l. 



(3.28) 



where is the permutation operator acting on the tensor product of the quantum 
Vj\ = Vj O Vj^ PI, Note that (13^2^ ) holds only for intermediate 



spaces as Fj^j[Vj,-^ 

sites with = j while the operators Rj-^^j and Rjj^^j corresponding to the end-points of the 
chain do not possess such property. As a consequence, the calculation of the Hamiltonian 
( p.27| ) deviates the standard derivation and requires a special consideration. 



^At this point we disagree with the statement made in iQ that the Hamiltonian ( 3.291 ) arises from 
the expansion of the fundamental transfer matrix. 
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Substituting ( |3.25| ) into ( p. 7] ) and ( |3.27| ) and making use of ( |3.28| ) we obtain after 
some algebra 



M~2 



k=2 



where the notation was introduced for the pair-wise Hamiltonian 



and the operator A7Y is defined as 

AH = tTj {Hj^j{uji) + Rj^j{-iuJi)Hj^j{0)Rj^j{iuJi) - Hj^j^{uji)} /tijl. 



(3.29) 



(3.30) 



(3.31) 



In comparison with the definition of the homogeneous spin chain Hamiltonian |T7], the 
expression ( |3.29| ) contains the additional operator A7i. One can show however that ATi 
provides a c- number correction to Ti^/. To this end one differentiates the both sides of 
the Yang-Baxter equation 



Rjji{u — iuJi)Rj2j{u)Rj2j-^{—iuJi 



R 



■nil ^ 



iuL)i)Rj^j{u)Rjj^{u — iuoi 



with respect to the spectral parameter u and puts m = 0. Then, multiplying the both 
sides of the relation from the right by and using ( p.28| ) one arrives at 

HjjA^i) + Rjj,{-iuJi)Hj^j{0)Rjj,{iiJi) = Hj^j,{u;i) + Rj^j^{uJi)Hj^j{0)Rj^j^{-iuJi) . 

(3.32) 

Using this identity one evaluates (|3.31| ) as0 



AH = RnU^u;i) ^''^'"i^^^ RnnH^i) = ^^^^^^ = const x 1 



tr,! 



32 



(3.33) 



In what follows we shall neglect the additive correction due to ATi. 

Thus, the integrable Hamiltonian of the inhomogeneous open spin chain with M sites 
is given by (|3.29|) with the pair-wise Hamiltonian defined by (|3.3CI|) and ( p.3|) as 



Hnj2{^) = ^(jl2 -uj)+ ^(ji2 +uj)- ipiJi + J2 - ^) - i'Ui + J2 + uj) 



(3.34) 



with iIj{x) = d\nr{x)/dx and the operator ji2 introduced in ( p.4| ). This expression is 
valid for arbitrary spins ji, j2 and, as it was anticipated, it is an even function of the 
impurity parameter u. 

It is now straightforward to see that for M = 3 and the parameters jk and Uk 
given by ( p.20| ), ( p.21|) and (|3.22| ), the Hamiltonian (|3.29| ) coincides (up to an overall 

^Here, the last relation in the r.h.s. follows from the SL{2) invariance of the 2-particle Hamiltonian, 
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normalization) with the QCD evolution kernels ( |2.1| ) and ( p.2| ). Following ( p.23| ) we 
introduce the interpolating Hamiltonian depending on the flow parameters cui^a 

n{uji,uj3) = i){ii2 - uji) + V^(ji2 + c^i) - - tui) - ^(| + t^l) 

+ V-fe - CU3) + ^(j23 + c^s) - ^(1 - CU3) - ^(1 + ^3) • (3.35) 

Then, the relation between 7i(co'i, cus) and the QCD evolution kernels, Eqs. ( p.l| ) and 
(PD, looks like 

Thus, the energies £3^ and entering into the solutions to the QCD evolution equations 
( p.l3| ) for the twist-3 quark-gluon distributions, Ds±{xi) and Dxixi), respectively, can 



be calculated as 

= ^{ll) + Y + ' = ^ ) + y + ' ^2-2^) 

where S{uji,lj3) are eigenvalues of the Hamiltonian ( p. 351 ) and l/A^"^— corrections corre- 
spond to nonplanar part of the evolution kernels (|2.1| ) and ( p.2| ). 

Having identified the QCD evolution kernels ( p.l|) and ( |2.2|) as Hamiltonians of the 
open spin chain model we now turn to the eigenvalue problem for ( p.29| ). The complete 
integrability of the Hamiltonian ( p.29| ) implies that TYa/ is a (complicated) function of the 
conserved charges. Using this dependence to which we shall refer as to the dispersion 
curve we can find the spectrum of the Hamiltonian by replacing the charges by their 
corresponding eigenvalues. However, instead of considering separately the eigenproblem 
for Hm we shall solve a more general problem of diagonalization the transfer matrices 
tj{u) of different spins j. 

We recall that the conserved charges and the Hamiltonian are generated by two 
different transfer matrices, t_i/2{u) and tj^{u), respectively. As we will show in Sect. 4, 
the transfer matrix t_i/2(w) can be diagonalized using the Algebraic Bethe Ansatz (ABA) 

. This immediately gives the eigenvalues of the integrals of motion. Moreover, as we 
will argue in Sect. 5, the transfer matrices of higher spins, tj{u), are related to t_i/2('u) 
through the nonlinear recurrence relations - the so-called fusion hierarchy |23|, ^ 



25| , Together with the ABA expressions for eigenvalues of t_i/2{u) these relations 
allow to reconstruct the spectrum of tj{u) and finally solve the eigenproblem for Hm- 

4. Diagonalization of the conserved charges 

Let us find the spectrum of the transfer matrix t_i/2{u). Throughout this section we 



shall keep the number of sites M arbitrary. According to the definition ( |3.17| ), t_i/2(M) 
is an operator acting on the quantum space of M particles, V,^ ® ... ® Vjj^j. The quantum 
space of the k—th particle, Vj^, has the highest weight 

\Ok) = l, L_,,jOfc)=0, LojjOk) = Jk\Ok) . (4.1) 
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Combining together the highest weights of M particles we construct the pseudovacuum 
state 



M 



(4.2) 



k=l 



with La being the total spin operator defined in (p.?]). As we will see in a moment, 
is the eigenstate of the Hamiltonian (|3.29|) with the spin = 0. 



4.1. Algebraic Bethe Ansatz 

The existence of the pseudovacuum state ( [4 .21 ) allows to apply the Algebraic Bethe 
Ansatz (ABA) for diagonalization of the transfer matrix t-i/2{u). Following the standard 
procedure [0, |1^, we first examine the action of the spin— 1/2 monodromy operator 
T_\/2{u) on the pseudovacuum state (^4.2|) . 

Due to the definition ( |3.5|) , T_i/2(m) is given by a product of the Lax operators (|3.1CI| ). 
Acting on the pseudovacuum state ( [4.2|) the Lax operators take the triangle form 



l.,Su-iUkm^)=('' + ^ . ° , .. )|^^+). (4.3) 

This allows to calculate their product in Eq. (|3.5|) as 

a{u)\n+) = 5ju+ , h{u)\n+) = , d{u)\n+) = 5(u+'-] , (4.4) 



where the operators a{u)...d{u) were defined in (|3.11| ) as different components of T_i/2(m) 
and the notation was introduced for the functions 

M 

6±{u) = \[{u-i{uJk±jk)). (4.5) 



k=l 



Substituting ( p. Ill) into (|3.16|) one gets 



T-i/2{u) = tr (^1^1 ^1^1 ) , U„{u) = A{u) + D{u) , (4.6) 
where the operators A{u) and D{u) are defined as 

Aiu) = a{ u ]d( —u I —biu |c( —u 

D{u) = d(^u - - - c(^n - - y (4.7) 

and the explicit form of the operators C{u) and B{u) is not relevant for our purposes. 
Using ( |4.4| ) together with ( |3.12| ) we find that the pseudovacuum diagonalizes both oper- 
ators 

A{u)\Q+) = 6+{u)6^{-u)\Q+) (4.8) 
D{u)\Q^) = |(^l- A^5^(-«)5_(«) + A5^(«)5_(-«)||f]+). (4.9) 
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and, as a consequence, we obtain from ( |4.6| ) 

U/2iu)\n^) = + ^5+(-«)<5-h} • (4.10) 

This relation implies that the pseudovacuum state is an eigenstate of the Hamiltonian 
( p.29| ). Comparing ([4 .21) with ( p.l4|) we identify as the eigenstate with the minimal 
conformal spin, = 0, 

^n=o{xi,...,xm) = \^+) = 1. (4.11) 

To find the corresponding eigenvalues of the conserved charges one has to expand the 
r.h.s. of (|4.1CI|) in powers of m^. In particular, for M = 3 and the parameters jk and 
Uk given by ( |3.20| ) and ( |3.23| ) we match (|4.10| ) into (|3.16|) to obtain the corresponding 
eigenvalue of the integral of motion as 

Qiuj„uj,)\n^) = f (f - ^1 - ^3) • (4.12) 

One can also get this relation using the definition (|3.23|) and taking into account that 
L'^f^\Q+) = {ji + jk){ji + jk — 1)1^+). The same identity allows to calculate the corre- 
sponding energy ( p. 351 ) as 

n{uJu uJs)'^N=o{zi, zm) = , (4.13) 

or equivalently S{lji,lj3) = at = 0. 

For higher values of the conformal spin > 1 the eigenstates of the auxiliary transfer 
matrix, or equivalently the eigenstates of the Hamiltonian Hm, are given by the Bethe 
States 113 

^;v(xi,...,xm) = |Ai,...,A;v) = C(Al)...C(A^)|^]+) (4.14) 

and are parameterized by the set of (complex) parameters {A^}. Here, C{u) is off- 
diagonal element of the monodromy operator defined in (f4.6|) . It becomes straightforward 
to verify that these states satisfy the highest weight condition ( p.l4| ) and diagonalize the 
transfer matrix provided that the complex parameters {Xj} satisfy certain conditions 
to be specified below in Eq. (|4.16| ). The corresponding eigenvalues of the normalized 
transfer matrix are given by 

2u jj^ {Xj+u){Xj-u) 

2u-i^ . X e- / X -TT (A7 + M + i)(A,- - M - i) /,,^\ 
2u [Aj + u){Aj — u) 

They differ from ( [4.10| ) by two "dressing" factors depending on the Bethe roots {Aj} 
and having poles at m = ±Aj. Since according to the definition ( p.l6| ), t_i/2(M) is an 
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even polynomial in u of degree 2M, the residue of the r.h.s. of ( [4.15| ) aX u = ±Aj and 



M = should vanish. One verifies that the pole at m = cancels in the sum of two terms 
( [4.15| ). The condition for t^i/2{u) to have zero residue aX u = \k leads to the set of the 
Bethe equations on the complex parameters {A^} 

(Afc + Aj+i)(Afc-Aj+i) ^ A+(Afc) 

where the notation was introduced for the functions 

A-,H = = (-1)*'^ - ^3t){u - 131) (4.17) 

fc=i 

A_(m) = A+(-m) 
with the parameters given by (|3.15| ). 

4.2. Baxter Q— function 

To calculate the spectrum of the transfer matrix ( |4.15|) one has to find all possible 
solutions Ai, Aa? to the Bethe Equation ( [4 .161) for different values of the conformal 



spin A^. Unfortunately, there exist no regular way of solving the Bethe Equations (|4.16| ). 
A more efficient (although equivalent) way of calculating the transfer matrix is based on 
the Baxter Q— function P7| . 

Let us define the Baxter Q— function as an even polynomial in the spectral parameter 
u of degree IN with the zeros given by the Bethe roots 

N 

Q(«) = n(" + A,)(«-A,). (4.18) 



Then, the obtained expression for the transfer matrix, (^4.15 ), is equivalent to the follow- 
ing 2nd order finite-difference equation on Q(u) 

t~\l2[u) Q(m) = A+(m) Q(m - 2) + A_(n) Q(m + i) (4.19) 

with A±('u) given by ( [4.17] ). Solving this equation one is interesting in finding all possible 
'polynomial solutions for Q(m). Each polynomial solution to the Baxter equation is in one- 
to-one correspondence with the Bethe states ( |4.14| ) - the degree of the polynomial defines 
the conformal spin of the states and its zeros provide the Bethe roots. In particular, the 
pseudovacuum state (|4.11| ) is related to a trivial solution to the Baxter equation 



Q;v=o(m) = 1. (4.20) 

We shall study the properties of the Baxter equation in more detail in Sects. 6 and 7. 

Thus, applying the Algebraic Bethe Ansatz we were able to diagonalize the transfer 
matrix t_i/2(M)- Expanding ( |4.15| ) in powers of u we can reconstruct the spectrum of the 
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integrals of motion. The next step should be diagonalization of the transfer matrix tj (u) 
generating the Hamiltonian ( |3.29|) . We recall that the only difference between t_if2{u) 
and t^j{u) is in the spin of the auxiliary space Vj entering the definition (|3.7|) . It is well 
known that higher (integer or half-integer) spin— j representations can be obtained from 
2j copies of spin— 1/2 representations through the so-called fusion procedure. Applying 
the same procedure to the transfer matrix of higher spin and performing fusion in the 
auxiliary space |jT8|] one can express t-j{u) in terms of spin— 1/2 transfer matrices t_ii2{u) 
and then generalize it to arbitrary j. Such relations - the so-called SU (2) fusion hierarchy 
- were first derived for the closed spin chains ||23| and later generalized to the open spin 
chains p4l I 
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5. Fusion hierarchy 

The fusion hierarchy provides the set of nonlinear recurrence relations between the trans- 
fer matrices of different spins. These relations have the same form for the closed and 
open spin chains pSfH 

f_j(n)Li/2(u + 2ij) = t_j-i/2{u) +t,j+i/2{u) f{u + (2j - l)z) , to{u) = 1 , (5.1) 

but the explicit form of the c- valued function f{u) is different in two cases. The fusion 
relation (|5.1|) is a direct consequence of the spin addition rule and the shifts of the 
spectral parameter in the arguments of t_i/2 and / are due to the fact that transfer 
matrices are built from noncommutative i?— operators. Here, t-j{u) is the normalized 
transfer matrix of the SU{2) spin j 

= t_j{-u - ij)pjiu) , (5.2) 
which differs from ( p. 71) by a multiplicative c-valued factor 

Pj{u) = Y\ p{u + ik) (5.3) 

k=0 



that compensates all poles of t_j{u). Similar to t_i/2(n), Eqs. (|3.13|) - (|3.17| ), the nor- 



malized transfer matrix t-j{u) is given by a polynomial in u of degree AMj with operator 
valued coefficients. 

To specify the function f{u) entering let us supplement the fusion relations ( |5.1| ) 
by the expression for t_i/2(w) in terms of the Baxter Q— function ( |4.19| ) 



r.V.M = A,M%lil + A_(,.)%±l). (5,4) 

Then, applying (|5.1|) and ( ^.4|) one can find the transfer matrices of higher spins through 
a single Q— function. In particular, for j = —1/2 one obtains from ( |5.1D the spin— 1 
transfer matrix as 

tll(M) = tli/2(M) t-i/2iu + i) - f{u) (5.5) 
^°Our definition of the SL{2) transfer matrix is related to that used in as t-j ~ T^^\ 
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and taking into account ( |5.4| ) one gets 



tli(n) = Q(n - i)Q{u + 2i) 



A_(m)A_(m 



Q(M-i)Q(M) Q(M)Q(M + i) 



+ 



A+(M)A+(M + i) 

Q(M + i)Q(M + 2i 

(5.6) 



with = A_(u)A+(u + — /(m). Then, the function f{u) can be fixed from the 

condition that t_i('u) should not have Q— independent term, 5t-i{u) = 0, leading to 



/(«) = A_HA_ 



(5.7) 



This expression is in agreement with the results of thorough analysis |2^, ^ . It is easy 
to see using (|5.1|) and (|5.4|) that the same condition ( ^.Tj) ensures that Q— independent 
terms do not appear in the expression for the transfer matrix of any spin.[^ Substitution 
of (|;T3) into d^) yields 

{2u - i){2u + 3i) 



fiu) 



Au(u + i) 



-5^{—u)5^{u)6+{u + i)6^{—u — i) . 



It is now straightforward to apply the fusion relations, Eqs.( |5.1| ) and ( |5.7| ), to calculate 
the transfer matrices of high spins. However, even though one could calculate in the way 
tj{u) for any given fixed j, it becomes difficult to find a general solution to the recurrence 
relations ( |5.1| ) in terms of t_i/2('u) which in turn are given by the ABA expression ( |4.15| ). 
Let us show that this problem can be avoided by using the Baxter Q— function. 



5.1. Eigenvalues of the transfer matrices 

Since the Baxter Q— function satisfies the second-order finite difference equation (|4.19| ), 
one expects to find its two linear independent solutions, Q+(m) and Q_(m). One of 
them, Q+(m), provides the polynomial solution ( 4.18 ) with the roots satisfying the Bethe 
equations. The second solution, Q_(m), is related to Q+(?/) through the Wronskian 
condition 

w{u) = Q+(m - i)Q-{u) - Q+(u)Q_(u - i) . (5.8) 
It follows from the Baxter equation ( [4.19| ) that w{u) satisfies the relations 

2j — l 

w{u + i) ^ A+{u) w{u + 2ij) ^ i-j- A+{u + ik) 

w{u) ~ A_(m) ' w{u) ~ j^J-^ A_{u + ik) ' 



As a trivial consequence of the definition (|5.8|) one gets the following useful identity 
Q_(M + 2ji) Q_(M-i) 



Q+{u + 2ji) Q+(m — i) 



(5.10) 



w{u) ^ w{u + i) ^ ^ w{u + 2ji) 



Q+(?i - z)Q+(n) Q+{u)Q+{u + i) Q+(m + (2j - l)z)Q+(u + 2ji) 



^^We would like to stress that this result is not sensitive to the explicit form of the functions A±(it) 
and therefore it holds for both open and closed spin chains. 
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Applying ( p.9| ) and ( [5.10| ) it is easy to verify that the spin— 1/2 and spin— 1 transfer 
matrices, Eqs. ( ^.41 ) and (|5^ ), respectively, can be expressed in terms of Q±(m) as 



t-l/2iu) 



t-Au) 



w{u) 



wiu 



[Q+{u - z)Q_(n + z) - Q+(m + i)Q_(M - i)] (5.11) 
[Q+(m - i)Q_(M + 2?) - Q+(m + 2i)Q_(M - t)] .(5.12) 



It now becomes straightforward to generalize these expressions to an arbitrary spin j 



w{u) 



[Q+{u - i)Q_(M + 2tj) - Q+(m + 2tj)Q_{u - t)] (5.13) 



which can be verified to be the exact solution to the fusion hierarchy (5J.). One can 
further simplify this expression by using (|5.1CI|) and writing it entirely in terms of the 
polynomial solution Q+('u) = Q(m) 



t.j{u) = Q{u-i)Q{u + 2ji) 



(5.14) 



X 



A(2i.o)| 



u 



+ ■ ■ ■ + 



_Q{u-i)Q{u) Q{u)Q{u + i) Q{u + {2j - l)i)Q{u + 2jt) 

Here the notation was introduced for the functions 



A(2i.o)(M) 
A(2^'-''')(m) 



2j-l 

JJ A_(m + ki) 

k=0 

2j-l l-l 

Yl A^{u + ki) Yl A+(m + 



rm] 



1 < ^ < 2 j - 1 



(5.15) 



k=l 



2j-l 



in=0 



(u + mi) 



m=0 



with A±{u) given by ( |4.17| ). 

Eq. ( |5.14D is our final expression for the normalized transfer matrix of the inhomo- 
geneous open spin chain. One verifies using ( [4.18| ) that for j = —1/2 it reproduces the 
ABA expression for the spin— 1/2 transfer matrix ([4.15|). 



5.2. Eigenvalues of the Hamiltonian 



Calculating the logarithmic derivative of t_j{u) at u = —ij and putting j = —j2 we 



can find the eigenvalues of the Hamiltonian by making use of Eqs. ( |5.2| ) and ( ^3.271 ). We 
notice from ( [4. 171 ) and ( |3.6| ) that A+(m) vanishes at u = —ij = ij2 leading to 



A+(-^j+e)=0(e 



2{M-2)^ 



(5.16) 
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As a consequence, one finds from ( [5.15| ) tliat A^^-' '''^ 



l,...,2j2 wliile A(2i'0)(-zj + e) = 0{e^) and tlie expansion of tlie normalized transfer 



matrix ( ^.141 ) around u = —ij = ij2 looks as 



Q(e - ij] 



A(2^'°)(-«j+e) + 0(e 



2(M-2)-i 



Finally, substituting this expansion into ( ^.21) and ( |3.27| ) and putting j 
late the energy of the open spin chain as 



S 



M 



Q'(ZJ2) Q'(-U2) 



Q(U2) Q(-«j2) 



(5.17) 
-j2 we calcu- 

(5.18) 



where j2 = ••• = Jm-i and prime denotes a derivative with respect to the spectral 
parameter. Here, the constant Em gets contribution from the p— factor in ( p. 21 ) as well as 
the constant term A7i in ( p.29| ) and determines the overall normalization of the energy. 
Since Em does not depend on the conformal spin of the state A^, it can be fixed using 
the expression for the energy of the pseudovacuum state ( [4.13|) . The latter is described 
by a trivial Baxter Q— function (|4.2CI| ). Substituting this solution into (|5.18| ) we obtain 
using ( CT ) 

Em = Sm{N = 0) = 0. 



Taking into account ( [1.18| ) we express the energy ( |5.19 ) as 



;WJ2 

' Q(«J2 



N 



2j2 



2 ' 



(5.19) 



where we indicated explicitly the dependence of the energy on the conformal spin and 



the number of sites M. Remarkably enough the relation ( |5.19| ) defining the energy of 
the open spin chain coincides with similar expressions for the energy of the closed spin 
chain ||19|, ^ although the Bethe roots and the Q— functions are different in two cases. 

Eqs. ( p.l4|) and ( p.l9|) are the main results of this section. We shall use them in the 
next section to find the spectrum of the Hamiltonian Hm- 



6. Baxter equation 

It follows from our analysis that the Baxter Q— function determines the spectrum of 
the transfer matrices, Eq. ( ^.141 ), and therefore plays a fundamental role in solving the 



eigenproblem for the open spin chains. Throughout this section we will be mostly inter- 
ested in solving the Baxter equation for the spin chain with M = 3 sites and parameters 
defined by the QCD evolution equations ( |3.2(]| ). However the methods described below 



are general enough and can be applied for solving the Baxter equation for arbitrary open 
spin chain. 
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6.1. Wilson polynomials 



Before analyzing the Baxter equation for M = 3 sites let us solve ([4.19 ) in the simplest 



case of the spin chain with M = 2 sites. This becomes useful for a number of reasons. 
Firstly, the M = 2 Baxter equation can be solved exactly. As we will show, its solutions 
can be identified as Wilson orthogonal polynomials p8[ and we shall use them as a basis 
for expanding the Baxter equation solutions for M > 3. Secondly, for the open spin 
chain with M = 2 sites the Hamiltonian ( |3.29| ) coincides with the pair-wise Hamiltonian 



( ^.34| ) and its diagonalization becomes trivial. As an independent check of Eq. ( |5.19| ) one 
should be able to obtain the same expression for the energy using the exact solution to 
the M = 2 Baxter equation. 

To write the Baxter equation (|4.19| ) for the spin chain with M = 2 sites one needs 



the expression for the transfer matrix ^-1/2(1*) defined in ( p. 161) . Going through the 



calculation of ( |3.16| ) and ( pj.llD for M = 2, we obtain the following relation 



t-i/2M=2{u) = A+(m) + A+(-m) - + l)NiN - 1 + 2ji + 2j2) , (6.1) 

where ji and j2 are the spins of the particles, ju = ji + j2 + N is the conformal spin of 
the state and A+(m) is given by (^4.17|) for M = 2 as 

2u -\- i 

A+(m) = {u - ijt){u - iji){u - ijt){u - ij^) ■ (6.2) 



Eq.( |6.lD allows to rewrite the M = 2 Baxter equation ( |4.19| ) in the form 



N{N + 2n + 2j2 - l)Q(n) = [Q{u - i) - Q(«)] + ^^TT [^(^ + ^) " ' 

Aw^ + 1 4^^ + 1 

(6.3) 

which has a striking similarity to the definition of the Wilson polynomials (see Eq.( p.3| ) 
in Appendix B). This immediately leads to 

Qm=2(m) = WV(W^; Ji^, jr, J2^, J2") = Wn{,U^]31 + ^l^jl - ^1,32 + t^2, J2 - . (6.4) 



We conclude that the polynomial solutions to the M = 2 Baxter equation ( |4.19| ) are given 



by the Wilson orthogonal polynomials with the parameters fixed by spins and impurities 
of the model. This leads to the following property of the Bethe roots, Eq. ([4.18|) . Being 



zeros of orthogonal polynomials they are real, simple and interlaced for different values 
of the spin N. At large N their distribution on the real axis can be described by a 
continuous distribution density P^ . 

To obtain the energy Sm=2 we substitute ( |6.4| ) into ( [5.19| ), identify u = ui — 002 and 



put UJ2 = for simplicity. Then, using the identity (|B.7|) we get the expression 
£2{N) = ^(ji +J2 + N + U;)- + 32 + u) 
+ V^(ji + i2 + N -u)- + 32 - uj) , 



which reproduces the spectrum of two-particle (inhomogeneous) Hamiltonian ( ^.34| ) for 

3i2 = 3i + 32 + N. 
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6.2. Three-particle Baxter equation 



Let us turn to the Baxter equation for the open spin chain with M = 3 sites and 
choose the spins ji = js = 1, j2 = 3/2 to correspond to the QCD evolution equations, 
Eqs. (13.201) , (|3.21|) and (|3.22|) . In addition we put 0^2 = and let the flow parameters ui 
and UJ3 to vary according to ( p.24| ). We obtain from (|4.17|) 



AM 



u 



U 1 

2 



J]^ (m - i{l + Uk)) {u - i{l - LOk)) 



(6.5) 



fc=l,3 



Evaluating ( p.l7| ) one can write the transfer matrix as 



4A^(A^ + 6) [u^ + 



t.i/2iu) = A4u)+A4-u)+[u' + 



2q- 15(cJi +u;|) + 



225" 



4 
(6.6) 

where integer N defines the conformal spin of the state ( |2.14| ), and q = q{uji,uj3) stands 
for the eigenvalue of the conserved charge ( p.23|) . These relations allow to represent the 
Baxter equation ( [4.19|) in the form 



A+(u) [Q(m - z) - Q(n)] + A+i-u) [Q{u + i) - Q{u)] 



(6.7) 



u 



2q - 15{u;l + ul) + 



225 



Examining the asymptotics of the both sides of this relation at large u one finds that 



Q(w) 



u 



2N 



. Together with the symmetry of ( |6.7| ) under 



u 



-u this implies that Q(m) 



is given by an even polynomial in u of degree 2N 



Q(n) = n2^ + aiM2A^-2 + ... + a^ 



N 
k=l 



U 



AD 



(6.8) 



with Ofc being some coefficients and \k being the Bethe roots. Substitution of this ex- 
pansion into ( |6.7| ) leads to an overcomplete system of linear equations on the coefficients 
Ofc whose consistency conditions give rise to the quantization conditions on the charge 
g(ci;i, cjs). Their solutions provide the set of quantized g(ti;i,ti;3) and the corresponding 
Baxter Q— functions. Finally, one applies (|5.19|) to calculate the energy as 



£M=3iN,q) = i 



.Q'(3V2) 



N 



Q(3V2) 



X] \2 , 9 • 



(6.9) 



fc=i 



6.2.1. Exact solutions 

Although we can not find the general solution to the Baxter equation (|6.7|) , there exist 
two special values of the flow parameter, (uijU^) = (3/2, 1/2) and {ui,uj2.) = (3/2,3/2), 
for which an additional degeneracy occurs and the Q— function can be calculated exactly. 
Remarkably enough it is for these two values of ui and uj^ that one recovers the QCD 
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evolution kernels ( p.24| ). The corresponding exact anomalous dimensions were first found 



in [^0|, ^ using different technique. 
Exact solution for Qs+- 

For (cuijCUa) = (3/2,1/2) it follows from (|6^ ) that A+(±m) acquires the same factor 
{u^ + 1/4) as the r.h.s. of (|6.7|). This allows to simplify ( |6.7| ) as 



A+(m) [Q{u - i) - Q(m)] + A+(-m) [Q(m + Q{u)] 
(4^2 + l)iV(Ar + 6) - 2(g - )! Q(m) , 



where 



A+(«) = - 
and the notation was introduced 



1 / I 
u \ 2 



3 

u 1 

2 



5 

U 2 

2 



„(0) 

Is 



N{N + 6) 



75 



Putting M = i/2 in (|6.10 ) and using the symmetry property Q(— m) 
that the l.h.s. of (|6.1CI| ) vanishes leading to 



(g-gf)Q(V2) = 0. 



(6.10) 
(6.11) 

(6.12) 
Q('u) one finds 

(6.13) 



This relation has two solutions: either q = q^g \ or Q{i/2) = 0. 

In the first case, substituting q = qf^ into ( |6.10D we notice that the Q— function 
satisfies the relation that is identical to the Baxter equation ( |6.3|) for the spin chain with 
M = 2 particles of spins ji = 3/2, j2 = 2 and impurities Ui = 0, uj2 = 1/2. This leads to 



Q 



(0) 

s+ 



u 



Wm u 



J 3 3 3 5 
' 2' 2'2' 2 



(6.14) 



Note that the Bethe roots of (|6.14| ) are real. Substituting Q^|(u) into (|5.19| ) and using 



( [B.Tp we find the energy of the "exact" level as 



£f+{N) = i;{N + 4) + i;{N + 3) - ^(4) - V^(3) 



2V'(iV + 3) + 



10 



(6.15) 



The explicit expression for the corresponding wave function can be found in ([A.6|) . 

In the second case, for q ^ qf^ the Baxter equation ( |6.10| ) can not be solved exactly. 
Nevertheless the relation Q(i/2) = ensures that its possible solutions have a pair of 
pure imaginary Bethe roots A± = ±i/2 



Qsiu) 



(6.16) 
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with P/v-i(w^) being a polynomial of degree — 1 in -u^, while the remaining 2N — 2 
roots are real. It follows from (|6.9|) that the pair of the Bethe roots A± = ±z/2 provides 
a contribution to the energy that is independent on the spin and that is bigger than 
the contribution of a real root Aq > by an amount 

A£ = ^ , - ^ , > 3 _ 1 ^ 1 _ (6.17) 

Al + 9/4 Ag + 9/4 2 3 6 ^ ' 

Since this contribution is present for all levels except the exact one, we expect that the 
energy of the exact level (|6.15|) should be separated from the rest of the spectrum by 
a finite "mass gap". As we will see in Sect. 7, this is indeed the case and ( |6.17D is in 
agreement with the large A^ calculation of the mass gap Eq. ( [7.41| ). 

Exact solution for Qt- 

For cji = = 3/2 we find from ( |6.5[ ) that A+(m) vanishes at m = — i/2 as A+(-u) ~ 
(m + Therefore, expanding the both sides of ( |6.7| ) around u = —i/2 up to 0{{u + 

i/2)^)— order we neglect the terms in ( |6.7| ) containing A+(-u) and obtain the system of 
two linear equations 

Q(z/2) [N{N + 6) - gr - f ] = 8zQ'(^/2) (6.18) 

Q(V2) [N{N + 6) + gT + f ] = iQ'{i/2) [N{N + 6) - + f ] ■ 

It has two different solutions: either the charge takes one of the following two values 

4+) = iV(Ar + 2) - ^ , q^-^ = NiN + 10) + ^ , (6.19) 

or Q{i/2) = Q'ii/2) = 0. 

In the first case, one is able to find the exact solutions to the Baxter equation. 
Similar to ( |6.14|) , these solutions are expressed in terms of the two-particle Q— functions. 



However in contrast with ( |6.14|) they are given by the sum of two Wilson polynomials of 
degree N and A^ — 1, respectivel}|3 

Q?^M = W^iv(«^;y,y)-a^W^^-i(«^;y,y) , (6.20) 



where 



-"^^It is interesting to note that the Qm=2— functions entering ( |6.14 ) and ( 3.20 ) are hnearly related to 
each other as 



Wn (m^; f , |, |, I) = a^WM (u^; §, |, |, I) + bNWN-i (u^; 
with aN = iN + 6)/{2N + 6) and 6^ = -N{N + 2){N + 3)/2. 
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We observe that the exact solutions ( |6.20| ) and ( |6.14| ) are related to each other as 



(N + 2) Q, 



(+), 



u 



NQ 



(-), 



u 



4(A^ + 3) 
N + 6 



Q 



(0), 



u 



(6.22) 



Examining the zeros of the polynomials ( |6.20| ) one finds that Q^'* and Q^"* have a pair 
of pure imaginary mutually conjugated Bethe roots and the remaining 2(A^ — 1) roots 
are real. 

Substituting ( |fci.2(JD into ( [5. 191 ) and taking into account the identities ( |.B.6D and ( [B.7| ) 
we find the corresponding exact energy levels as 



S^-\n) = 2^{N + 3) + 



h 27, 

N + 3 3 

^ 8 



(6.23) 



Note that £^^^''(0) = in accordance with ( 4.13 ) and the level spacing, 



'(+)/ 



Sj, (iV) = —4/ (A^ + 3), vanishes at large A^. The property of the Q— functions (|6.22|) is 
translated into the following relation between the energies 



2 

3' 



(6.24) 



which is valid for arbitrary spin A^. 

As we will see in Sect. 7, the exact solutions, ( |6.15| ) and (|6.23| ), are the lowest energy 
levels in the spectrum of the corresponding Hamiltonians. In addition, invariance of the 
Hamiltonian Ht = 7i(3/2, 3/2) and the charge Qt = (5(3/2, 3/2) under permutations of 
the end-points allows to assign a definite parity to the energy levels. One verifies using 
the explicit expression for the wave function, Eq. (|A.6| ), that two exact energy levels, 



s!^^ and 



(-) 



.rp diiu i^rp have an opposite parity under xi 



X3 that alternates as A^ changes. 



The second solution to the system (|6.18| ), Q,(i/2) = Q'(i/2) = 0, implies that the 
Q— functions describing "nonexact" levels, q ^ \ have a double degenerate pair of 
pure imaginary Bethe roots ±i/2 



u 



4 I Pn-2{u'' 



(6.25) 



Similar to ( |6.16| ), this property leads to an appearance of the gap separating the exact 
levels ( |6.23|) of the Hamiltonian Tix from the rest of the spectrum. 



6.2.2. Master recurrence relations 

To find the rest of the spectrum of the Hamiltonian 7i(ci;i, cus) we shall use the Wilson 
polynomials as basis for solving the Baxter equation ( |6.7|) . Apart from the fact that 
this basis is very convenient due to its orthogonality and completeness properties, the 
expansion of the M = 3 Baxter Q— functions over the M = 2 solutions has a simple 
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interpretation. It is in one-to-one correspondence with the decomposition of the 3-body 
wave function over the 2-particle states. 

More exphcitly, the wave function of the M = 3 particle state with given conformal 
spin h = ji + j2 + js + N can be decomposed over the set of states having the same 
total spin h and, in addition, possessing a definite conformal spin in the channel defined 
by any pair of particles. For instance, the decomposition of the eigenstate ^'tv in the 
(12)— channel looks as 

TV 

^7v(x.) = 5^c„<f (x,), (6.26) 

n=0 

where ^I/^vn^ the conformal spin ji2 = ji + j2 + n in the channel (12) 

^L^Sf (^0=Jl2(jl2-l)v&Sf (X.) 

f 12)3 

and c„ are the expansion coefficients. The states ^l/\rn automatically diagonalize the 
two-particle Hamiltonian, ?ii2 and one associates with them the corresponding Wilson 
polynomials Qn^'*^(u) = Wniu"^)- Then, the decomposition of the M = 3 solution of the 
Baxter equation takes the same form as ( |6.26| ) 



N 



Q^(«) = 5^c„Qf)^(tz) (6.27) 



ra=0 



with c„ being the same coefficients as in (|6.26|) . The reason for this is that (|6.26|) and 



( |6.27| ) describe the same eigenstate but in different representations. The Q— function 



determines the wave function in the separated coordinates and it is related to ( p.26| ) 
through a unitary transformation well known as the Separation of Variables ||32|| . 

1 ( 23) 

Choosing another pair of particles we find different polynomials Qn yu) that are 

f 12)3 

linear related to Qn (u) through the Racah 6j— symbols of the SL{2) group. We obtain 
their explicit expressions by replacing the parameters ji and uji entering ( |6.4| ) by their 
values corresponding to the (12)— and (23)— pairs of particles with ji = ja = 1, j2 = 3/2 
and a;2 = 0. In this way we get 

Q(^^)3(n) = Wn {u'; H 1 + ^1, 1 - ^i) (6.28) 
3)(^) = H 1 + 1 - . (6.29) 



Substitution of (|6.27|) into ( |6.7D leads to the system of recurrence relations on the 



expansion coefficients c^. To find their explicit form it proves convenient to define the 
following functions 

A+(«) = -^±^^, Q(«)=(^2 + .;|)Q(«). (6.30) 
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Here the additional factor was introduced to remove the contribution to ( |6.5| ) of the 
particle with the impurity tus. It brings A+(n) to the form corresponding to the spin 
chain with M = 2 sites, Eq. ( |6.2|) , and the Baxter equation (|6.?1 ) can be rewritten as 



4^2 + 1 



A, 



-u 



4m2 + 1 
{N+l){N + 5)Q{u 



Q{u + i) — Qiu] 

q-qo{N) 



with 



^ 3 — Aiu'n 35 o 3 o 
qo{N) = N{N + 6)^-^ - -ujj - -ujj + 



105 



Q(m) (6.31) 



(6.32) 



Note that the l.h.s. of this equation coincides with the two-particle Baxter equation 
( |6.3[ ) in the (12)— channel and the total two-particle spin ji2 = ji + j2 + N + 1 = | + -/V. 
According to ( |6.4| ), the solutions to this reduced equation are given by basis function 
Qiv+i (■") defined in (|6.28|) . Considering the r.h.s. of (|6.31| ) as a perturbation we seek for 
a general solution to ( |6.31| ) in the form 



N+l 



Q(«) = E^"Q 



(12)3, 



U 



(6.33) 



k=0 



Substituting this expansion into (|6.31| ) and taking into account the Baxter equation for 
Qr^'(M), Eq. (U, we get 

N N+l 



Y,{N+l-k){N + 5 + k)m 



(12)3, 



U 



(12)3, 



U 



(6.34) 



k=0 



k=0 



Note that according to ( |6.33|) and (|6.3CI|) the r.h.s. of this relation is given by {u^ + 
ul)Q{u). This leads to the following expression for Q{u) (up to an irrelevant factor) 



N 



QH = + l-k){N + 5 + k)fkQ 



(12)3, 



U 



(6.35) 



k=0 



which coincides with ( |6.27| ) after an appropriate redefinition of the expansion coefficients. 
Applying the identity (|B.4|) and using orthogonality of the Wilson polynomials we find 
from (|6.34|) that the coefficients fk satisfy the three-term recurrence relations. To simplify 
their form it is convenient to replace fk by a new set of coefficients Uk defined as 



Uk = fk{N + l-k){N + 5 + k) Qf ^'(3^/2) 



with )'(3V2) 



(3)fc (I + a;i)fc (I 



uji)k due to Eqs. 



(6.36) 

and ([B.6|). Then, the 



expansion ( |6.35| ) takes the form 

N o(12)3 



Q 



u 



N 



k=0 



(3)fc(| + cJi)fc(| -cji)fc 



k=0 



2 + iu, 2 
3, 9 + 1^1, 9 — 1^1 



-k,k + 4,^ + iu, \ — iu 



1 

(6.37) 
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with the coefficients Uk satisfying the three-term recurrence relations 

-Uk = Uk-iAk-i - Uk \ Ak + Ck + ujr^ - - ] + Uk+iCk+i 



2{N + l-k){N + 5 + k) . y • . • ,i 

(6.38) 

Here qo{N) was defined in ( |6.32D and the coefficients Ak and Ck are given by Eqs.( [B.5|) 
with a = b = 3/2, c = 1 + ui and d = 1 — ui 

{k + 3){k + A){k + I + uji){k + I - cui) 



A, 



Ck 



{2k + A){2k + 5) 
k{k + l){k + I + uji){k + \-uJi] 



{2k + ?>){2k + A) 

Substituting ( |6.37|) into ( |6.9|) and using the identities ( p.6|) and (p.Tp we find that 
the expression for the energy takes a remarkably simple form 

X _ Ef=0 ^fc(g) [^(fc + i + ^l) + V'(fc + I - Ul) - ^(1 + Ul) - ^(1 - g^l)] 

(6.39) 

where we indicated explicitly the dependence of the energy on the spin and the 
conserved charge q. This expression has a simple interpretation. The total energy is given 
by the sum over all possible conformal spins k in the two-particle (12)— channel. The 
contribution of the spin k is proportional to the two-particle energy SM=2{k) weighted 
with the factor Uk/ J2n '^n- 

The master recurrence relations ( |6.38| ) allow to construct the solutions to the Baxter 
equation ( |6.37| ) and, as a consequence, find the spectrum of the conserved charge q and 
the energy S{N,q). Solving ( |6.38| ) we impose the normalization condition UQ^q) = 1. 
Under this choice the coefficients Un{q) are given by polynomials in q of degree n. For 
n = -|- 1 one finds from (|6.38|) that 

UN+i{q)=0. (6.40) 

This relation provides the quantization condition on the charge q. Solving ( |6.40D we 
get -|- 1 different solutions for q that have the properties of roots of the orthogonal 
polynomial: quantized q are real, nondegenerate and interlaced for different values of 
N. Then, it follows from ( |6.38|) and (|6.39|) that the expansion coefficients Uk{q) and the 
energy £{N, q) are also real. 

The results of the numerical solution of the recurrence relations for the charge q and 
the energy £{N, q) and their dependence on the spin and the impurity parameters cji 3 
are shown in Figs. |I| and |^, respectively. We observe from Figs. ^ that the eigenvalues 
of the charge q and the energy £{N, q) grow with the spin A^ and occupy the band 

gmin(A^) <q< gmax(A^) , ^mm(A^) < S{N, q) < S^^^N) (6.41) 

with gmin = 0{N^), S^in ~ 21nA^ and gmax = O^N"^), £^max ~ 41nA^. In addition, the 
spectrum of q and S{N, q) exhibit an interesting properties of regularity. To understand 
these properties we shall apply in the next section the asymptotic methods to find an 
approximate solution to the Baxter equation for large values of the spin A^. 
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Figure 1: The dependence of the quantized values of the charge q and the energy £(N,q) on 
the spin N for the values of the impurity parameters uJi =003 = 3/2 corresponding to the 
evolution kernel Tir- The dotted lines represent two trajectories with i = (upper curve) and 
£ = 3 (lower curve) described in Sect. 7.4. 
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6.2.3. Singular states 



According to (|6.38|) tlie solutions to the recurrence relations depend on the shift param- 
eters uJk- Examining the dependence of the energy ( |6.39|) on uji and uj^ shown in Fig. ^a 



one observes that S{N,q) diverges as ook approaches half-integer values 

The total number of the singular points is equal to the spin A^. As A^ increases they 
start to occupy the whole real axis except of the interval 

- ^ < ^1,3 < I , (6.43) 

which defines the stability region of the model. 

This phenomenon can be understood by noticing that the two-particle energy in 
the (12)— channel entering ( |6.39| ) diverges as Uk casing due to singularities of the 
V'— functions. Taking ui = 5/2 + n with n = , 1 , ... we find from ( |6.39| ) that divergence 



comes from the contribution of two particle spins A; = n -|- 1, A^. As a consequence, 
for this value of Ui the only levels in the spectrum that have a finite energy are those 
with Uk = for A; = n + 1, A^. If this condition is not satisfied the energy becomes 
infinite.0 For ui = 5/2 + n the total number of finite energy states is equal to n + 1 
and at 77, = A^ there are no divergent states left in the spectrum. Since the finite energy 
states form an orthogonal subspace, the singular states have Uk/u^ = for = 0, n. 

According to ( |6.9|) , the existence of the singular states corresponds to appearance of 
the Bethe roots at = ±3'i/2. 



7. Asymptotic solution of the Baxter equation 



Although the Baxter equation ( |6.7| ) can not be solved exactly there is a simple way to 
find its asymptotic solutions at large values of the spin A^ The method is based 

on the observation that the both sides of the Baxter equation ( |4.19| ) have a different 
scaling behaviour at large A^. We find that for u = fixed the l.h.s. of ( [4.19|) scales as 
A±{u) ~ whereas its r.h.s. involves the transfer matrix (|6.6| ) that grows at large A^ 
as 

U/2{u) = 0{q) (7.1) 
depending to the value of the charge q{N), Eq. ( |6.41| ). Then, introducing the function 



Qiu) 



(7.2) 



13 



Even though singularities cancel in the numerator of ( |6.39 ) under a weaker condition ^ 



+ 1 "fe 



examining the singular part of the matrix element of 



0, one finds a stronger condition X)fc=n 
two-particle Hamiltonian in the (12)— channel, {'^ n\Hi2\'^ n) with 'Sat given by ( 6.26| ) 
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Figure 2: The flow of the energy £{N, q) and the conserved charge q = 2N{N + 6)(z^^ + 3/4) 
with 1V3 for N = 5 and uJi = 3/2. Two vertical dotted hnes correspond to 0^3 = 1/2 and 
UJ3 = 3/2 and define the spectrum of the evolution kernels Tis and TCt, respectively. 



one rewrites the Baxter equation ( [4. 191) 



+ A+(-m) (p{u) = t-i/2{u) . 



(7.3) 



Lp{u — i) 

A general solution to this relation is given by (an infinite) continuous fraction. Taking 



into account ( |7.1| ) we find that at large N and u = fixed the solutions to ( |7.3| ) are given 
by 



-1/2 



u] 



AJ-u] 



+ 0{l/q). 



(7.4) 
(7.5) 
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It is easy to verify that these relations lead to the following expressions for the Q— function 



Q-(w) = Qh 



r(iu) 



T{i{u + ak))T{i{u - ak)) 



n 



-u 



(7.6) 
(7.7) 



where were defined in ( ^.ISj ) and (7^ denote the roots of the transfer matrix ( p.6[ ) 



-1/2 



-2ni 



M — (Tt 



(7.8) 



k=l 



The general solution to the Baxter equation is given by a linear combination of Q±(m). 
Requiring Q(m) to be an even function of u one gets 



QUu) = Q+{u) + Q-{u) = Q+(m) + Q+(-m) . 



(7.9) 



We would like to stress that thus defined asymptotic Q— function obeys the Baxter 
equation ( |6.7| ) only in the limit N —>■ oo and u = fixed up to (9(l/g)— corrections. 

The dependence of the Q— function on the charge q enters into ( [7.9|) through the 
roots cTfc of the transfer matrix, Eq. (|7.8|) . Since large scales, q and A^(A^ + 6), appear in 
( |b.6D with a common prefactor (m^ + 1/4), one of the roots is g— independent and it is 
given by 

al = -] + 0{N-^) . (7.10) 



Then, taking into account the values of the parameters one simplifies ( |7.6| ) as 



Q+iu) 



,»r(m)r(zu+ 2; 



n 



T(iu - iak)Tiiu + iak] 



r^(^M + I) r(m + 1 - uJk)T{iu +l + uJk) 



(7.11) 



fc=l,3 



The values of the remaining roots (Xi^s vary significantly as the charge q changes inside 
the band (ICTD . 



Close to the upper bound q = 0{N^) one finds that the both roots are real and 
increase with the spin as 



N{N + 6) + 0{N^) 



(7.12) 



<3 = [1 ± VI - ^ 

where q = q/{N{N + Q)f = 0{N^). 

Close to the lower bound q = 0{N'^) one finds that one of the roots is A^— independent 

(tI = 2N{N + Q) + 0{N^) , al = v"^ + 0{N-^) , (7.13) 

where the parameter v is defined as 



V 



2N{N + 6) 4 



- = 0{N" 



(7.14) 
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Notice that v takes real values for q > 3N{N + 6)/2 and it becomes pure imaginary for 
q < 3N{N + 6)/2. This happens in particular for the exact levels, Eqs. (|6.12| ) and ( p.l9| ) 



(4±))2 = -\t^ + 0{N~^) , 4 = -l + 0{N-^) . (7.15) 

As we will see later, the appearance of the complex roots of the transfer matrix is closely 
related to the existence of the mass gap in the spectrum of the Hamiltonian. 



7.1. Dispersion curve 



Let us apply the solution (|7i9|) to obtain the asymptotic expression for the energy 8. 
According to ( |6.9| ) the energy is determined by a logarithmic derivative of the Q— function 
at M = ±3i/2. These values of u belong to the applicability region of the asymptotic 
solutions, N ^ oo and u = fixed, and therefore one is allowed to replace Q— function in 
( |6.9|) by its asymptotic expression ( [7.9| ). 

It follows from ( [7.11|) that Q-|-(±3z/2) = and close to u = ±3z/2 the Q— function 
is given by one of the functions, Q = Q±(u). Substitution of (|7.9|) into (|0|) yields 



(7.16) 



£UN,q) = [^(l + ^^fe) -In2-£(cui,cu3) + 0(l/g) 

fe=l,3 

where the notation was introduced for the normalization constant 

e{uji, ujs) = ^ ^ (I + ^k) + ^ (I - ^k) 

fc=l,3 



(7.17) 



and (Ti_3 are nontrivial roots of the transfer matrix defined in Eqs. (|7.8| ) and ( |6.6|) . The ex- 
plicit expressions for ai^s are quite cumbersome and one can use instead the approximate 
expression 



a 



1,3 



1 ± Vl-g N{N + 6) + 0(A^°) , q 



q-lN{N + 6) 
[N{N + 6)]2 



(7.18) 



which agrees with the asymptotic expressions ( |7.12| ) and ( [7.131 ). The relation ( [7.1(j| ) 
establishes the dependence of the energy on the integrals of motion, N and q, and 
provides an explicit expression for the dispersion curve of the spin chain model. 

The expression ( |7.16 ) simplifies in the upper part of the spectrum ( |6.41 ). We find 
using (|7l^) and ( ^l3D that 



SUN, q) = ln(g/2) - 6{ui, U3) + 0{N- 

provided that q = 0{N^). 

In the lower part of the spectrum we get from ( |7.13 ) 



(7.19) 



SUN, q) = HNiN + 6)]+^{ + i;{--tu] - 5(0^1, cus) + 0{N-^) (7.20) 
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with = q/{2N{N + 6)) — 3/4 = 0{N^). For the special values of the impurity pa- 
rameters, Eqs. ( |3.21| ) and ( p.22| ), the relations ( [7.19| ) and ( [7 .201 ) coincide with analogous 
expressions obtained in [l^. 

Since u could be either real or pure imaginary in ( [7. 201) , we separate the corresponding 



energy levels into two groups by introducing the "vacuum" energy 

S^UN) = HN{N + 6)] + 2^ (f ) - e{uju oos) . (7.21) 

The energy levels with z/^ > lie above the vacuum, Sas{N,q) > S^g^dN), and we 
shall refer to them as "continuum", while a few lowest energy levels with z/^ < and 
£as{N,q) < £^vac(^) correspond to the "bound states". 

As we will show in Sect. 7.2.1 (see Eq. ( [7.30| ) below), the parameters scale in the 
continuum at large as z/ ~ 1/ In A^. For small real u one expands the r.h.s. of ( [7.20| ) 
in powers of u and finds that the energy in the continuum grows linearly with q close to 

^vac 

q 3" 



^continuum {N,q)-S,UN) 



7C(3)) . (7.22) 



_A^(A^ + 6) 2_ 

The energy of the bound states has a completely different behaviour at large A^ because 
in contrast with the continuum the parameter |z/| takes finite values for these states. As 
a consequence, their energy levels are separated from the continuum by a gap 

AS = S^UN) - ^bound(iV, q) = 2^ (0 - ^ Q + kl) - ^ Q - l^^l) • (7-23) 

We will estimate the value of the mass gap at large A^ in Sect. 7.2.3. 

The asymptotic expansions ( 7.16 ), ( 7.19 ) and ( [7. 201 ) are valid up to corrections van- 



ishing at large A^. Let us check the accuracy of (|7.16|) by comparing Sas{N,q) with 
the exact values of the energy obtained from the solution of the recurrence relations, 
Eqs. ( |6.39|) and ( |6.3^ ). To this end we choose the shift parameters to correspond to 
the QCD evolution kernel ( p.22|) , ui = us = 3/2, and take the spin to be A^ = 10. 



Applying (|6.38|) and ( |6.39|) we find A^ + 1 pairs of the exact eigenvalues (g. Sex) that we 
compare with the dispersion curve Si^{N,q) as shown in Fig. ^. Calculating the differ- 
ence 6S = Sa_s{N,q) — Sax{N,q) we find that the asymptotic formula ( |7.16| ) reproduces 
the exact energy with a very high accuracy 6S < 10~^, which increases up to 6S < 10~^ 
as one goes from the lower part of the spectrum to the upper bound by increasing q. 

7.2. Quantization of the conserved charge 

Let us show that the matching of the analytical properties of the asymptotic solutions 
( [7.9| ) into Eq. ( p. 81 ) leads to the quantization conditions on the charge q. 

It follows from ( |7.9| ) and ( [7.11| ) that the function Qas(w) has an infinite series of 



poles in the complex m— plane generated by F— functions in the numerator of ( [7.11 ) and 
located at 

Mpoie=|±^, n, ±ai + in, ±as + iri^ (7.24) 
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q/ (N(N+6) )^ 

Figure 3: The dependence on the energy £{N,q) on the conserved charge q at N = 10. 
Crosses denote the exact values of the energy and the sohd Hue corresponds to the asymptotic 
expression ( [7.16|) . 



with n being arbitrary integer. This seems to be in contradiction with the fact that the 
exact solutions to the Baxter equation are given by polynomials in u, Eq. ( |6.8|) . Notice 
however that the asymptotic and the exact solutions for Q(m) should coincide (up to 
corrections vanishing at large A^) only in a finite region of u whereas some poles in ( [r.24| ) 
are moved outside this region as — *• cxo. This allows to remove from the consideration 
the "moving" poles in ( [7.24| ) satisfying the condition Upoie = 0{N) and keep only the 
"fixed" poles, Wpoie = C(A^°). In particular, in the upper part of the spectrum we get 
from (|7.12|) that ai^s = 0{N) and therefore the fixed poles of Q(m) are located at 



''^fixed pole 



in} . (7.25) 



Similarly in the lower part of the spectrum we use ( [7.13|) and ([7.241) to find the fixed 
poles of ( [7.11| ) as 



^fixcd pole 



q=0(N'2) [ 2 

with 1^ defined in ( [7.14 ) and n = integer. 



±- , in, ±i' + in\- (7.26) 



Let us now examine the residue of Qas(^) at the fixed poles (|7.25|) in the upper part 
of the spectrum. Using the explicit expression for the functions Q+('u), Eq. ( [7.11| ), one 
finds that at large the residue of Q+('u) at u = im with m = |,1,2,... is suppressed 
with respect to its residue at m = by a factor (cricia)""^™' = 0{q~'^"^) that vanishes as 
N oo. Therefore up to 0(l/g)— corrections one can neglect in ( [7.25] ) the fixed poles 
at M = ±1 and u = in and keep only the leading pole at u = 0. 

Repeating similar consideration for the fixed poles in the lower part of the spectrum, 
( [7.26| ), one finds that only the poles at m = and u = ±z/ survive in the large N limit 
provided that |lmz/| < 1/2. Notice that the residue of Q(m) at m = ±z/ scales as ~ g^Si;^ 
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and for |Imi/| > 1/2 it becomes comparable with the contribution of subleading poles at 
u = ±i/2 and u = in ^ 0.0 Thus, in order to match the asymptotic solution ( [7 .91) into 
( |678D one has to require that Q(m) should have a zero residue at u = 0. In addition, in 
the lower part of the spectrum one has to impose the same condition at m = 

It is easy to see from (|7.6| ) that the residue of Q+(m) and Q-('u) at m = cancel each 
other in the sum ( |7.9| ) independently on the value of the charge q. At the same time, 
calculating the residue of ( [7 .91 ) and ( |7.11| ) at u = ±zz/ we find after some algebra that it 
vanishes provided that u obeys the following equation 

where v'^ = q/{2N{N + 6)) — 3/4. This relation establishes the quantization condition 
on the charge q in the lower part of the spectrum q = 0{N'^). 

Solving the quantization conditions ( [7.27| ) we distinguish two cases: the states in the 
continuum, u"^ > 0, and the bound states, < 0. 

7.2.1. Solving the quantization conditions in continuum 

For > one gets from (|7.27] ) 



V In NiN + 6) - arg ^^'=^'' \, , ^ = ni (7.28) 

with i being an integer. For the special values of the impurity parameters, Eqs. ( ^.21| ) 
and ( |3.22| ), this relation coincides with the quantization conditions obtained in [^, 



It follows from ( [7.28|) that u scales as ~ c/ ln[iV(A^+6)] at large with the constant 
c having a nontrivial dependence on the impurity parameters ui and u^. In particular, 
varying these parameters one finds that c gets a finite contribution as cui^a passes through 
integer values. To show this we use the identity 

arg r(l — cUfc + iu) = arg r(2 — Uk + w) — arg {1 — uJk + w) (7.29) 

and notice that as uj^ passes through the value tUjt = 1 the phase of (1 — cjjt — iv) changes 
by TT. This transition occurs in the region |1 — cj^l = 0{v) = 0{1/ \n[N{N + 6)]). In 
a similar manner, changing Uk from 1 to oo we find that every time Uk passes integer 
positive values the l.h.s. of ( [7.28 ) increases by vr, or equivalently the integer i in the r.h.s. 



of ( [7.28|) decreases as — 1. As we will see in a moment, this effect corresponds 

to formation of the bound state which "dives" below the vacuum energy and causes 
reparameterization of the levels in the continuum. 

For < ciJi , < 1 and 1 — cui^a ^ 1/ ln[A^(A^ + 6] we obtain the solution to ( |7.28| ) as 

ln[N{N + 6)] + C ^ [h^N, 

^^One can systematically improve the accuracy of the asymptotic solution by including nonleading 
1/q corrections to the asymptotic solutions ([/.(j). 



^oW = TZT^TTTT-l^^r-T^ + ( ) (7-30) 
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with 

C = 4^(1) - 2^(1) - [^(1 - ^k) + ^(1 + ujk)] . 

k=l,3 

Substituting ( [fISOD into (|O0|) we find the energy of the £— th level as 



f,^£.. + (14C(3)-16);-,i^^|^l^ + 0(^) (7.31) 

with the vacuum energy S^^c given by ( [7.21D . Integer i = 1,2,... enumerates the energy 
levels with i = 1 being the lowest level in the continuum. Note that this expression is 
valid only for the lowest / ~ ln[A^(A^ + 6)] states. Using ( |7.31| ) we find the level spacing 
in the continuum close to the vacuum level as 

S£dN) O ( / . (7.32) 



Let us now increase uji and keep ^3 unchanged. We find from (|7.31|) that the energy 



of all levels increases (see Fig. ||a) while the difference S£ — E^^c decreases with ui (see 
Fig.^). The lowest level £ii=i rapidly approaches the vacuum energy ^vac as cJi 1. 
Once uii crosses the value tui = 1 one applies ( |7.29| ) to obtain the solution to ( [7. 281) as 

K^)L<.,<. =^0(^-1) (7.33) 

0<U3<1 



with £ > 1 and t'o(^) given by (|7.3CI| ). Thus, the lowest energy level (i = 1 disappears from 



the continuum and this induces reparameterization of the remaining levels, — 
Since this effect occurs every time as uji or uj^ passes through integer values it is now 
easy to write the general solution to ( [7.28| ) valid for arbitrary values of uji^^ 

= z/o(^-N-N), (7.34) 

where [uJk\ denotes an integer part of uok- We conclude that as a result of the flow of the 
spectrum of the model from cji = = to arbitrary u;i^3, 

A/Lound = [t^i] + [t^s] (7.35) 

states cross the vacuum level with i/vac = and the energy £^vac and disappear from the 
continuum. These levels have the charge z/^ < 0, or equivalently q < 3N{N + 6)/2, and 
they can be described using ( [7.27| ). 



7.2.2. Solving the quantization conditions for the bound states 

For i/^ < we put u = ip with p > and examine the large behavior of the both 
sides of (|7.27| ). The l.h.s. of (|7.27|) grows as ~ N'^^ whereas the r.h.s. is a meromorphic 
A^— independent function of p. Therefore in the limit N ^ 00 one could satisfy ( [7.271 ) 



38 



either through a trivial solution p ~ 1/ ln[A^(A^ + 6)], or taking p to be close to the poles 
of the r.h.s. 

We find using (|7.27| ) that for p > and < uji^z < 1 the poles are located at p = 1/2 



and p = positive integer. However, these poles coincide with the subleading fixed poles in 
( [7.26| ) and their contribution is beyond an approximation at which (|7.27] ) was obtained.]^ 



Increasing the values of the impurity parameter cui > 1 we observe that the r.h.s. of 
( [7.27] ) develops the poles at 

p = cui — 1 , tui — 2 , ... , {uji — 1} (7.36) 

with {uoi — 1} being a fractional part of uoi. Similar phenomenon occurs as one increases 
the parameter uo^. Then, for arbitrary positive uoi and uj^ the total number of the poles 
is equal to [cji] + [cus] and coincides with the number of the "missing" levels from the 
continuum Abound- Thus, at large N Eq. ( [7.27] ) has two branches of the solutions each 
parameterized by integers rii and ns, respectively 



nu + {o^fc - 1} + O (iV-i) (7.37) 



with = 1, 3 and ra^ = 1 , ... , [tUfc— 1]. To find (9(1/A^) corrections to these expressions one 
has to include subleading corrections to the quantization condition ( [7.27|) , or equivalently 
to the asymptotic solution ( [7. 91 ) and ( |7.ti| ). Each solution ( [7. 37] ) corresponds to the bound 
state with the energy £^bound('^i) or £^bound('^3) given by 

i^vac -inbound (^fe) = 2^ Q + ^fc + {^k " 1}^ -'^ ~ ~ i^'^) 

(7.38) 

with < rik < [cUfc — 1]. Thus, as uJk passes through an integer the level from continuum 
( [7.31| ) with the energy £^£=„j. crosses the vacuum to transform into the bound state with 



the energy inbound ('^fc)- This transition occurs in the region \u!k — nk\ = 0(1/ In A^) and 
it is clearly seen on Fig. |a. For given Ui and cua the total number of bound states is 
equal to ( [7.351 ). Note that up to 0{1/N) corrections the difference £^vac — ^^bound does not 
depend on the spin N. 

In particular, using ([7.35[) and ( [7.37[) we find that for ui = 1/2 , = 3/2 and ui = 



UJ3 = 3/2 the spectrum contains A/bound = 1 and Abound = 2 bound states, respectively, 
with iu = —1/2 + 0{1/N). In the second case, two bound states are degenerate up to 
0{1/N) corrections. This is in agreement with the exact results ([6.23|) and ( |6.15[) . 



7.2.3. Mass gap 

We observe from Figs. |l|b and ^jb that at large and fixed 001^3 the bound states are 
separated from the continuum by a finite mass gap, A{lji,lj3) = min {Si — inbound)- We 
estimate its value using ( [7.31[) and ([7.38[ ) as 

A(u;i, cug) = £vac " max„^,„2 i^bound("'fc) + 0{1/ In^ A^) . (7.39) 

^^It follows from the matching condition that subleading corrections to the asymptotic Baxter equation 
solution should screen all poles in ([7.26). 
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Figure 4: The difference £{N,q) — <fvac(-^) for the lowest five energy levels of the Hamilto- 
nian TC{uji,uj3): (a) the dependence on the impurity CJ3 for N = 100 and cJi = 3/2; (b) the 
dependence on the spin for ui = ^3 = 3/2. 



Then, A(a;i,a;3) = for < ^^13 < 1 and 

A(^l, ^3) = 27/> (^j -^l;f^ + {uj,- 1} 



^ ( - - {^1 - 1} 



(7.40) 



both for < Co'3 < 1 < cui and for 1 < cui < UJ3. 

Applying (|7.40|) for the values of parameters to'1^3 corresponding to the QCD evolution 
kernel, (|3.22|) and (p.21|) , one calculates the mass gap as 



A 0, = A 0, = 2^ (I) - ^(1) - ^(2) = 3 - 41n2 = 0.227411 . (7.41) 

This result agrees with the exact calculations of the difference £{N,q) — Sva.c{N) at 
cji = ti;3 = 3/2 as shown in Fig. ^. Note that ( [7.41| ) is valid up to corrections ~ 
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-k"^ /\y?[N{N + 6)e'-^] which decrease slowly with and provide a sizeable contribution 
to A at finite A^. 

According to (|7.41|) the mass gaps in the spectrum of the QCD evolution kernels, T^t 



and 7^5+, are the same. Nevertheless, the energies of bound state for these two Hamil- 
tonians are different due to dependence of the vacuum energy on the shift parameters, 
Eqs. ( |7.21| ) and ( |7.17| ). Their difference is given by 



^•T,bound ~ ^S,bound — ^vac ( 2' 2 / ~ ^^'^ V 2' 2 / ~ 3 (7.42) 



and it coincides with the exact result ( |6.24 ). 



To summarize, matching the asymptotic solution ( |7.9| ) into the exact form ( |6.8[ ) at 
large N and u = fixed we obtained the quantization conditions ( |7.27| ) and used them to 
calculate the asymptotic expansion of the charge q and the energy S{N,q) in the lower 
part of the spectrum, Eqs. (|7.31|) and ( [7.38| ), and estimated the value of the mass gap. 



Eqs. ( [7.40|) and ( [7.41|) . To find similar expansion in the upper part of the spectrum 



one needs an asymptotic solution to the Baxter equation that is valid at large N and 
u = 0{N). 



7.3. WKB expansion 

Analyzing the Baxter equation ( |4.19| ) at large and u 
introduce the scaling variables 



0{N) it is convenient to 



u 



r]x, r]=N{N + 6) 
Then, one gets from ( |675|) and (|6.6| ) 



q = q/{N{N + 6))' 



with s = 1/2 -2^^jk = -13/2 and 

t_y2{vx) = 2r]^x^ [Vo{x) - 1 + 0{ri' 



Vo = 2/x^ - q/x^ 



(7.43) 



(7.44) 



(7.45) 



Substituting these relations into (|4.19 ) one observes that in the leading rj —>■ oo limit the 
Baxter equation takes the form of the discretized one- dimensional Schrodinger equation 
on the "wave function" Q(a;?7). In this equation the parameter l/rj plays the role of the 
Planck constant and Vo{x) defines the potential. This suggests to look for its solutions 
in the WKB form Q 

Q+ixT]) = exp {ir]S{x)) , S{x) = So{x) + r]-'Si{x) + C(r/-2) , (7.46) 

where So{x), Si{x), ... are r^— independent and the expansion is assumed to be uniformly 
convergent. Substituting the WKB ansatz (|7.46| ) into the Baxter equation ([4.19|) one 
expands its both sides in powers of l/r] to get 



cosh 5*0(0;) 



Voix), S[ix) 



-Sq{x) coth5*o(a;) — i— 



(7.47) 
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This leads to the following WKB expansion 



Q+{xri) 



A/sinh Sq{x] 



with Sq{x) satisfying ( |7.47| ). Finally, combining together ( [7.9| ) and ( |7.48| ) we arrive at 
the asymptotic expansion of the Baxter equation solution that is valid at large N and 
u = xr] = 0{N). 

Performing the matching of (|7.9|) and (|7.48|) into the exact solution (|6.8|) , we shall im- 
pose the following two conditions: the asymptotic solution Qwkb {xt]) should be a single 
valued function of x and it should oscillate on the real x— axis in order for Q— function 
to have the real Bethe roots = 0{N). 

Exploring the interpretation of Qwkb (a;//) as a quantum- mechanical WKB wave func- 
tion we obtain in the standard way that the first requirement is translated into the 
Borh-Sommerfeld quantization condition on the eikonal phase 



dx S'(x) 



(7.49) 



where / is a nonnegative integer and the integration is performed over a closed contour 
on the complex x— plane encircling the classical interval of motion X- < x < x^ defined 

as 

(1 - Vo{x±)f = (1 - 2x^2 + qxl'^Y = 1 . (7.50) 
Replacing S = Sq + Si/r] and using ( [7.47|) we get from (|7.49|) 



dx Sq{x) 



(7.51) 



Since the WKB wave function oscillates on the interval we can satisfy the second 

requirement by demanding x_ and x+ to be real. This leads to the following constraint 
on the possible values of the charge q = q/[N{N + 6)]^ 



< g < 



1 



(7.52) 



Moreover, solving the Borh-Sommerfeld quantization conditions ( [7.51|) one can develop 
the asymptotic expansion of the charge q close to the upper bound, q = 0{N^) 



q = q/[N{N + 6)f 



(1) 
V 



v(2) 



+ 



7]^ 



+ 



g(3) 

^3 



o(4) 



(7.53) 



with q^''^ being the expansion coefficients depending on integer I and the shift parameters 
a;i^3. Taking the limit 77 — *■ cxo in ( [7.51| ) one finds that the l.h.s. of (|7.51| ) has to vanish. 
This happens when the classically allowed region shrinks into a point, x+ = x_, or 
equivalently 

= i (7.54) 
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due to ( fT^Ol) . 

To calculate the first nonleading coefficient g^^-* we have to find preasymptotic 0{r]~^) 
term in the large r^— expansion of the l.h.s. of ( [7.51|) and match it against 0{ri~^) in the 
r.h.s. of (|7.51| ). In this way one gets 



(1) 



(7,55) 



Note that the coefficients ( [7.541) and (|7.55|) do not depend on the shift parameters a;i^3. 

The WKB expansion (|7.48|) can be systematically improved by taking into account 
additional terms in ( [7.4(j| ) and using the Baxter equation to express them in terms of 
Sq{x). Going through this procedure and applying the technique developed in 
were able to extend the expansion ( |7.53|) up to 0{ri~^) term 



we 



,(2) 



,(3) 
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The asymptotic expansion (|7.53| ) describes the spectrum of the conserved charge q close 
to the upper bound q = 0{N^). It follows from ( [7. 531 ) that quantized values of the 
charge form the trajectories q = q{i, N) labeled by a nonnegative integer i. An example 
of such a trajectory for £ = is shown in Fig. |I|. For the special values of the impurity 
parameters, Eqs. ( |3.21|) and ( |3.22|) , the first three terms of the expansion (|7.53|) agree 



with analogous expressions obtained in |14 

For the trajectories with i = fixed the expansion in ( [7. 531 ) goes over inverse powers 
of 1] N while for trajectories with i = 0{N) the coefficients q^^'^ grow as ~ A^*^. In the 
latter case, the WKB expansion can be improved by introducing a new parameter 

1 



+ 



V 



-V2. 



Then, reexpanding q in powers of l/rj for ^ = fixed one finds 
with 
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Here, in contrast with ( [7.53| ) the expansion goes only over even powers of t] and it 
describes the levels with I = 0{N). Note that the leading term go(0 does not depend 
on the impurity parameters and it can be found exactly as a solution to the Whitham 



equations 



The obtained WKB expansions ( [7.53|) and ( [7.571 ) combined with the asymptotic ex- 
pression for the dispersion curve ( |7.19| ) provide a good description of the spectrum close 
to the upper bound. 

7.4. Trajectories 

The relations (|7.30|) , (|7.14|) and (|7.53| ) define two different sets of the trajectories pa- 
rameterized by integers C. and £, respectively. Being combined together they provide a 
complimentary description of the spectrum throughout the continuum. Positive integer 
^ enumerates the trajectories ( [7.31| ) lying above the vacuum level in the lower part of 
the spectrum, while integer £ = 0,1,2,... corresponds to the ordering of the levels from 
above in the upper part of the spectrum. Since for fixed spin N the total number of the 
levels is equal to + 1, the integers ^ and £ are formally related to each other as 

I=N-l. (7.58) 

The examples of I— and £— trajectories for the Hamiltonian Tix are shown in Fig. |I[ 
The trajectories with £ = and £ = 3 go through the states in the continuum with 
the maximal amd minimal energy, respectively. Two trajectories lying below the £ = 3 
trajectory correspond to the bound states. 

One finds from ( 7.53| ) and ( [7.55 ) that the distance between two neighboring trajec- 



tories in the upper part of the spectrum behaves at large as 

5q=[N{N + Q)f5iq = 0{N^). (7.59) 

This expression should be compared with the level spacing in the lower part of the 
continuum that one finds using ( [7.14| ) and ( [r.30| ) as 

5q = 2N{N + Q)5i u^ = 0{ ^-^^ j . (7.60) 



Using the properties of the spectral curve, Eq. ( [7.19| ) and ( |7.31[ ), one obtains the corre- 
sponding energy level spacings in the continuum as 



5£{N, g) O ( ^) , 6£iN^ g) o(^). (7.61) 



We recall that the bound states are separated from the continuum by a finite mass gap 
i VM ) and iflM) . 



It is interesting to observe that the same relations ( [7.61[ ) describe the spectrum of 
the anomalous dimensions of the baryonic distribution amplitudes [^]. This suggests 
that ( [7.61| ) are universal features of the three-particle evolution equations in multi-color 



44 



QCD. We refer to for further discussion of properties of the trajectories and their 
physical interpretation. 

Introducing the trajectories one can classify the conformal operators (|2.1CI| ) for dif- 
ferent q as belonging to the different trajectories. Each trajectory describes a sepa- 
rate component of the twist-3 quark-gluon distribution (|2.13|) . In contrast with the 
/X— dependence of the distribution, the scale dependence of its components is of the 
DGLAP-type, Eq. ( |2.13| ), with the anomalous dimensions given by ( |2.11| ). Their mixing 
with other components is protected by the additional Q— symmetry of the model. 

In general, the quark-gluon distributions D{xi; fi) enter into physical observable in- 
tegrated over parton momentum fractions (see Eq. ( |A.7|) ). As a consequence, it gets 
contribution from all trajectories and it scale dependence becomes nontrivial. However 
depending on the particular form of the corresponding weight function one may encounter 
the situation when most of the trajectories decouple and only one trajectory contributes 
in the multi-color limit. In this case, the scale dependence significantly simplifies and 
takes the standard DGLAP form. As shown in the Appendix A, this is exactly what 
happens for the twist-3 nucleon structure functions in the multi-color limit. We would 
like to note that this property holds only in the leading iVc — > oo limit and the spin struc- 
ture functions get contribution from all trajectory through nonplanar l/A*"^— corrections. 
These corrections will modify the form of the twist-3 quark-gluon states constructed in 
this paper but will not destroy the analytical properties of the trajectories. 



8. Conclusions 

In this paper we studied the evolution equations for the twist-3 quark-gluon parton 
distributions in the multi-color QCD. The evolution equations follow from the scaling 
dependence of nonlocal twist-3 quark-gluon string operator on the light-cone and have 
the form of one-dimensional Schrodinger equation for three particles on the light-front. 
Our analysis was based on the observation that in the multi-color limit the evolution 
equations possess an additional integral of motion and turned out to be effectively equiv- 
alent to the Schrodinger equation for integrable open Heisenberg spin chain model. The 
parameters of the model are uniquely fixed by properties of underlying quark-gluon 
system. Using this correspondence we constructed the basis of local conformal twist-3 
quark-gluon operators and calculated their anomalous dimensions as the energy levels of 
the open spin magnets. We identified the integral of motion of the spin chain as a new 
quantum number that separates different components of the twist-3 parton distributions. 
Each component evolves independently and its scaling dependence is governed by the 
anomalous dimensions of the conformal operators. 

To find the spectrum of the QCD induced open Heisenberg spin magnet we developed 
the Bethe Ansatz technique based on the Baxter equation. Solving a nonlinear fusion 
relations for the transfer matrices of the (inhomogeneous) open spin chain models we 
derived the exact expression for the energy, or equivalently the anomalous dimensions 
of quark-gluon distributions, in terms of the Baxter Q— function. The properties of 
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the Baxter equation were studied in detail and its solutions were constructed using 
different asymptotic methods. We demonstrated that the obtained solutions provide 
a good qualitative description of the spectrum of the model and reveal a number of 
interesting properties: the fine structure of the energy spectrum is described by the set 
of trajectories, few lowest energy levels are separated from the rest of the spectrum by 
a finite mass gap, for certain values of the impurity parameters the energy diverges and 
the system becomes unstable. 

We believe that the open spin chain models define a new universality class for differ- 
ent problems in high-energy QCD like the Regge asymptotics of quark-gluon scattering 
amplitudes and the evolution equations for high-twist quark-gluon distributions and the 
results obtained in this paper could be applied there as well. 

Our consideration was restricted to the multi-color limit Nc oo. Nonleading 
corrections destroy integrability of the QCD evolution kernels, Eqs. ( |2.1| ) and 
( p.2| ), and modify the spectrum of the anomalous dimensions. Numerical calculations 
indicate 0, ^ that l/N^ corrections do not destroy the analytical structure of the 
trajectories but modify the level spacing in the upper part of spectrum generating a 
mass gap separating the highest energy level from the rest of the spectrum. These 
effects can be systematically taking into account following the approach |^ and deserve 
additional studies. 

Acknowledgments 

We are most grateful to V.M. Braun for illuminating discussions and for collaboration 
at the early stage. One of us (G.K.) would like to thank V.P. Spiridonov for useful 
discussions on the Wilson polynomials. This work was supported in part by DFG Project 
N KI-623/1-2 (S.D.) and by the EU network "Training and Mobility of Researchers", 
contract FMRX-CT98-0194 (G.K.). 

A Appendix: Twist-3 nucleon parton distributions 

In this appendix we summarize the relations between twist-3 nucleon structure functions 
and quark-gluon distributions introduced in (|TTT]). 

Following p we define the chiral-odd, e{x; fi) and h^^x^fi), and the chiral-even, 
gT{x]iJ,), structure functions through the matrix elements of nonlocal light-cone quark 
operators 

{p,s\q{y)q{-y)\p,s) = 2M j dx e^'^^^^y^ e{x] ^i) 

(p,s|g(y)yW75g(-2/)b,s) = 2{sy)M^ j\x^^^^y^hL{x;ii) (A.l) 
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{P, s\q{yhal5q{-y)\p, s) = 2Ms^a j dxe^'^'^Py'^grix; , 

where = zn^ with = and non-Abehan phase factor is omitted for brevity. Here, 
\p,s) is the nucleon state with momentum p^, p^ = and the spin s^, = —1. In 
addition, s±a denotes two-dimensional transverse component of the spin orthogonal to 
plane defined by the vectors and p^. 

Using the QCD equations of motion one can separate the twist-2 Wandzura-Wilczek 
contribution to the structure functions hL^x) and grix) |TU|, |^, ^ ^ 



dy 

hiix) = 2x I -^hi{y) + hiix) 

y 



X 



9t{x) = [ ^gM+M^) (A.2) 

Jx y 

with hi{y) and gi{y) being twist-2 nucleon distributions and and grix) being 

genuine twist-3 part of the structure functions. They can be expressed as pOl, pT| 



'^dxe'"^P^^e{x;fi) = ^j\u j\t{p,s\T^{u,t,-u)\p,s) 

dxe^^"(P"^/^^(x;/x) = -I ^^^l [ duu [ dtt{p, s\Ti^,{u,t, -u)\p, s) (A.3) 
1 [sz)M Jq J _^ 

X (P, s\{u - t)S^{u, t, -u) + {u + t)S~{u, t, -u) \p, s) . 

The matrix element entering the r.h.s. of these relations define (auxiliary) twist-3 quark- 
gluon distributions [§] 

(p,s|T|(M,t,-M)|p,s) = -AM{pnf j r'a;e'(^'")["^("-*)-"«("+*)]De(xi,X2,a;3) 
{p,s\Ti^,{u,t,-u)\p,s) = 2iM\sn){pn) [ e'(P")["i("-*)-"^("+*)]D,, Jxi, xa, X3OA.4) 



{p,s\S^{u,t,-u)\p,s) = 4:S±aM{pnf J I)a;e'(P")["i("-*)-"^("+*)lD±(xi, X2, X3) . 

Here, Vx = dxidx2dxz5{xi + X2 + X3) with xi, X2 and (— X3) being the momentum 
fractions carried by quark, gluon and antiquark, respectively. 

Expanding the both sides of (|A.3|) and ( |A.4|) in powers of {pn) we find the moments 
of the structure functions as 



dxx^ e{x) = / VxDeixij'^^^^^i^ii^^) 
J -I 
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^dxx^hLix) = j^-^ J ^ 
1 _ 1 /-i 

,(±,0) 



(A.5) 



L'+(Xi)^S5-2(^l'^3) - D_{Xi)'^'}^'_2i^3,Xl 



(0) 



where ^^jy ' ( homogeneous polynomials in xi and X3 of degree iV defined as 



'^^n\xi,X3) 



X 



N+1 



-X3 



Xi + X3 

X^ 



3:3 



Xi + X3 
N+2 



(A.6) 



Xi +X3 

Inverting these relations and taking into account the spectral properties of the distribu- 
tions one gets [0, [1^, ||, |3^ 



1 dx' 
e{x) = - I -Ye{x,x') 



f -\ Ju Jj 



dx' dx" 
hdx) = -x I — ——— [d,, - d,.] Yh{x',x") 



1 ^^1 /•! J™'/ 



(A.7) 



9t{x) 



r dx' r dx ' 

/ -r -^[d.'Y,ix',x") + d^.Y,ix",x')] 

J X X J _i X X 



where the notation was introduced for (anti)symmetrized twist-3 quark-gluon distribu- 
tions 



x" —x' 



F/,(x', x") = Dh{x', -x' + x", -x") - Dh{x", x' - x", -x') 

Yg{x',x") = Z)+(x', -x' + x", -x") + D„(x",x' - x", -x') . 

In contrast with the Dg.h.g— functions these distributions are real functions of the parton 
fractions. 

One can verify by a direct calculation that the polynomials ( |A.6| ) diagonalize the 
conserved charges Qt and Qs+, Eqs. ( |2.17| ) and ( |2.16| ), on the subspace Xi + X2 -l- X3 = 



(±)| 



Xj 



{±),T,(^)| 



Xi 



(0),T,(0)| 



Xi 



(A.8) 



with the eigenvalues q!^^ and g^^"* given by the exact solutions, Eqs. ( |6.19| ) and (|6.12| ), 
respectively. Then, comparing ( |A.5| ) with ( ^.13| ) we conclude that in the multi-color limit, 
Nc 00, the moments of the nucleon structure functions, Aif{N; jj) = J^^ dxx^ f{x; /i), 
have a scale dependence of the DGLAP type 



^^^Mf{N■, /i) = -^j{N)Mf{N; /i) + 0{1/N^) 



f = e,hL,gT 



(A.9) 
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with the anomalous dimensions equal to the energies of the exact levels, Eqs. ( |3.36| ) 
( |6.15| ) and (|6.23| ). For the chiral-odd structure functions one gets [^| 



asNc 
asNc 

277 



4-^)(iV-2) + ^ 



13 



4-)(Ar-2) + ^ 



13 



+ 0{l/Nl 



(A.IO) 



and for chiral-even structure function one finds the flavor nonsinglet (NS) contribution 
to the anomalous dimension as ||3C 



ifAN) 



27r 



£f\N~2) + 



II' 



+ 0{l/Nl 



(A.ll) 



This simplification occurs because among all possible components of the twist-3 quark- 
gluon distributions the twist-3 nucleon structure functions "select" only those corre- 
sponding to the lowest exact levels. The scale dependence of the chiral-odd structure 
functions e(x) and h^lx) is associated with two "exact" trajectories of the Hamiltonian 
T-Ct, while the chiral-even structure function grix) is governed by the "exact" trajectory 
of the Hamiltonian Tis- 



B Appendix: Wilson polynomials 



In this Appendix we summarize the properties of the Wilson polynomials p8[ . They are 
defined in terms of hypergeometric series as 



TT7" / '2 7 1\ / 1 \ / \ / 1\ 7~t / ^ ^ -^5 ^ tX ^ CI %CC 

Ty„(a; ;a,6,c,d) = (a + 6)„(a + c)„(a + d)„ 4F3 , , , , , 

a + o,a + c,a + d 



fB.ll 



with s = a + b + c + d, (a)„ = r(a + n)/r(a), n integer and a...d being arbitrary 
parameters. Wn are given by polynomials in x"^ of degree n and the prefactor in (|B.1|) is 



chosen to ensure the symmetry with respect to any permutations of the parameters 

Wnix"^; a, b, c, d) = Wn{x'^; b, a, c,d) = ... = Wn{x'^; a, b, d, c) . (B.2) 

The Wilson polynomials form the system orthogonal polynomials on a line —00 < x < 00 
and the corresponding orthogonality and completeness conditions can be found in . 
They satisfy the second-order finite difference equation 

n(n + s — l)y{x) = B{x) [y{x + i) — y{x)] + B{—x) [y{x — i) — y{x)] (B.3) 

with y{x) = Wn{x^] a, b, c, d) and 

^. . {x + ia){x + ib){x + ic){x + id) 
~ 2x{2x + i) ■ 
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In addition, Wn obeys the following three term recurrence relations 

X^Pnix^) = Pn+l{x') + {An + - a^)Pn{x') + A^-lCnPn-lix') • (B.4) 

Here, 

Pn{x') = i~ir ^'~l}'' W4x'-a,b,c,d) 



S - l)2n 



are normalized polynomials and 



^ {n + s — l){n + a + b){n + a + c){n + a + d) 
~ {2n + s- l)(2n + s) 

^ _ n{n + b + c-l){n + b + d-l){n + c + d-l) 



Using (|B.1|) one can obtain the following useful relations 

Wnix'^; a,b, c, d) = (a + 6)„(a + c)„(a + (B.6) 



x=±ia 

„2. 



±idx\jiWn{x ■,a,a,c,d) = ^jj{n + a + c) — 'ilj{a + c) (B.7) 

+ '?/'(n + a + c?) — ilj{a + c?) . 
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